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Abstract. Integrable Killing tensors are used to classify orthogonal coordi- 
nates in which the classical Hamilton-Jacobi equation can be solved by a sepa- 
ration of variables. We completely solve the Nijenhuis integrability conditions 
for Killing tensors on the sphere and give a set of isometry invariants for the 
integrability of a Killing tensor. We describe explicitly the space of solutions 
as well as its quotient under isometrics as projective varieties and interpret 
their algebro-geometric properties in terms of Killing tensors. Furthermore, 
we identify all Stackel systems in these varieties. This allows us to recover the 
known list of separation coordinates on in a simple and purely algebraic 
way. In particular, we prove that their moduli space is homeomorphic to the 
associahedron K^. 
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1.1. The problem. Given a partial differential equation, it is a natural problem 
to seek for coordinate systems in which this equation can be solved by a separation 
of variables and to classify all such coordinate systems. It is not surprising, that 
for classical equations this problem has a very long history and is marked by the 
contributions of a number of prominent mathematicians. For example, consider the 
Hamilton-Jacobi equation 
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on a Riemannian manifold. We say that it separates (additively) in a system of 
coordinates xi, . . . , x„ if it admits a solution of the form 



where each function Si only depends on the variable Xi. If this is the case, the 
initial partial differential equation decouples into n ordinary differential equations. 
We call the corresponding coordinate system separation coordinates. 

It turns out that the Hamilton- Jacobi equation in the form (1.1) separates in a 
given orthogonal coordinate system if and only if it also separates for V = and if 
the potential V satisfies a certain compatibility condition. In fact, the Hamilton- 
Jacobi equation with V = plays a key role for separation of variables in general. 
The reason is that for orthogonal coordinates the (additive) separation of this equa- 
tion is a necessary condition for the (multiplicative) separation of other classical 
equations such as the Laplace equation, the Helmholtz equation or the Schrodinger 
equation. Sufficient conditions can be given in the form of compatibility conditions 
for the curvature tensor and for the potential. 

The classification of orthogonal separation coordinates for the Hamilton-Jacobi 
equation with V = thus leads to a variety of ordinary differential equations. 
For Euclidean space, these are well-known equations and define special functions 
appearing all over in mathematics and physics - such as Bessel functions, Legendre 
polynomials, spherical harmonics or Mathieu functions, to name just a few. In 
particular, these functions serve as bases for series expansions in the explicit solution 
of boundary value problems on domains bounded by coordinate hypersurfaces. 

Separation coordinates are intimately related to the existence of rank two Killing 
tensors that satisfy a certain integrability condition. 

Definition 1.1. A symmetric tensor Ka/j on a Riemannian manifold is a Killing 
tensor if and only if 



where V is the Levi-Civita connection of the metric g and the round parenthesis 
denote complete symmetrisation of the enclosed indices. Depending on the con- 
text, we will consider a Killing tensor either as a symmetric bilinear form Kafj or 
alternatively as a symmetric endomorphism if"^ = g^^K^p. 

We say a symmetric endomorphism K is integrable, if around any point in some 
open and dense set we find local coordinates such that the coordinate vectors are 
eigenvectors of K. This geometric definition can be cast into an analytic form with 
the aid of the Nijenhuis torsion of K, defined by 
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Integrability in the above sense is then equivalent to the following system of nonli- 
near first order partial differential equations [Nij51]. We take this characterisation 
as the definition of integrability. 
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Definition 1.2. A symmetric endomorpliism field K"'^ on a Riemannian manifold 
is integrahle if and only if it satisfies the Nijenhuis integrability conditions 

O^N\p^g^]s (1.4a) 
= A^V^"1^ (1.4b) 
= N\^^K^^,K's, (1.4c) 

where the square brackets denote complete antisymmetrisation in the enclosed in- 
dices. 

At the end of the 19th century, Paul Stackel showed that any system of or- 
thogonal separation coordinates gives rise to n — 1 linearly independent, mutually 
commuting integrable Killing tensors with simple eigenvalues [Sta97]. And in 1934, 
Luther P. Eisenhart proved the converse [Eis34] . Notice that the metric is a trivial 
example of an integrable Killing tensor and it commutes with any other Killing 
tensor, but does not have simple eigenvalues. 

Definition 1.3. A Stackel system on an n-dimensional Riemannian manifold is an 
n-dimensional space of mutually commuting integrable Killing tensors. 

It can be shown that any integrable Killing tensor with simple eigenvalues is 
contained in a Stackel system and vice versa, any Stackel system contains such 
a Killing tensor [Ben97]. This establishes a one-to-one correspondence between 
orthogonal separation coordinates and Stackel systems. 

The Nijenhuis integrability conditions (1.4) are homogeneous algebraic equa- 
tions in K and V K, which are invariant under isometrics. Separation of variables 
therefore naturally leads to the following problem. 

Problem I. Determine the projective variety K.{M) of all integrable Killing tensors 
on a given Riemannian manifold M . That is, solve explicitely the Nijenhuis inte- 
grability conditions (1.4) for Killing tensors (1.2). Moreover, find natural isometry 
invariants characterising the integrability of a Killing tensor. 

By the definition of integrability, any Killing tensor in dimension one and two 
is integrable. Therefore, Problem I is trivial in dimensions less than three, once 
the space of Killing tensors is known. However, it has never been solved in higher 
dimensions, not even for the least non-trivial examples, namely Euclidean space 
M.^, the sphere or hyperbolic space H^. Moreover, an explicit solution of the 
equations (1.4) has been considered intractable by various experts in this domain 
[HMS05]. 

Although a general solution of Problem I is unknown for any non-trivial exam- 
ple, there is a construction which yields a large family of Stackel systems. This 
family is interesting in itself, because it appears under different names in a number 
of different guises in the theory of integrable systems: as "Newtonian systems of 
quasi-Lagrangian type" [RWML99], as "systems admitting special conformal Kil- 
ling tensors" [CSTOO, CraOSa], as "cofactor systems" [LunOl, RW09], as "bi-quasi- 
Hamiltonian systems" [CS02, CraOSb], as "L-systems" [Bcn05] or as "Benenti sys- 
tems" [BM03]. We follow the latter nomenclature of Bolsinov and Matveev, because 
their description is best suited to our context. This brings us to the following ques- 
tion. 

Question I. Do Benenti systems parametrise the entire variety of integrable Killing 
tensors ? 
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Moreover, by construction, the set of Benenti systems is not only invariant under 
isometries, but even under the projective group. So, if the answer to Question I 
is "no" and the isonietry group is a proper subgroup of the projective group, the 
following question arises. 

Question II. Is the variety of integrable Killing tensors invariant under the pro- 
jective group? 

Obviously, there are three groups acting on the variety /C(M) of integrable Killing 
tensors: 

(i) M, acting by addition of multiples of the metric g. 

(ii) The isometry group Isom(M). 

(iii) M*, acting by multiplication. 

The actions (i) and (iii) result from the fact that adding the identity to an endo- 
morphism or multiplying it by a constant does not alter the eigenspaces and hence 
does not affect integrability. The actions (i) and (ii) commute and the action (iii) 
descends to the quotient. Hence the essential information about integrable Killing 
tensors is encoded in the projectivisation of the quotient 



This leads to the following problem. 

Problem II. Determine explicitely the quotient (1.5), which describes integrable 
Killing tensors modulo obvious symmetries. 

Recall that orthogonal separation coordinates are in bijective correspondence 
with Stackel systems. Thus, in view of the classification problem for orthogonal 
separation coordinates, one is interested in linear subspaces of commuting integrable 
Killing tensors, rather than in the integrable Killing tensors themselves. This leads 
to the following problem. 

Problem III. Determine explicitely the space of all Stackel systems on a given 
Riemannian manifold M (modulo isometries). That is, find all maximal linear 
subspaces of commuting Killing tensors in K,{M). 

The separation coordinates can be recovered from an integrable Killing tensor 
with simple eigenvalues, because the coordinate hypersurfaces are orthogonal to 
the eigenvectors. This leads to the final step in the determination of separation 
coordinates via integrable Killing tensors and Stackel systems. 

Problem IV. Find a (canonical) representative in each Stackel system for which it 
is possible to solve the eigenvalue problem and which has simple eigenvalues. Deduce 
the classification of orthogonal separation coordinates and give natural isometry 
invariants. 

In this article we will completely solve Problems I to IV for the three dimen- 
sional sphere equipped with the round metric and give answers to the questions 
I and II. This is done without relying on any computer algebra. One of the reasons 
we choose is that it has the simplest isometry group among the 3-dimensional 
constant curvature manifolds. Our solution will reveal diverse facets of the problem 
of separation of variables, related to representation theory, algebraic geometry, geo- 
metric invariant theory, geodesic equivalence, Stasheff polytopes or moduli spaces 
of stable algebraic curves. 




(1.5) 
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1.2. Prior results. In his fundamental paper Eisenhart also derived a list of or- 
thogonal separation coordinates for M.^ and [Eis34]. These results have been 
generalised by Kalnins & Miller, who gave an iterative diagrammatical procedure 
for the construction of orthogonal separation coordinates on spaces of constant 
curvature in any dimension [Kal86, KM86]. Furthermore. Eisenhart computed the 
Stackel systems on M.^ in their respective separation coordinates. Despite these 
advances, the topology or the geometry of the space of integrable Killing tensors or 
its quotient by isometrics has never been understood for any non-trivial example. 

Since explicit solutions were out of sight, the use of computer algebra increased in 
recent years. This has lead to the construction of a set of 15 isometry invariants for 
Killing tensors on M.^ which discriminate between the different types of orthogonal 
separation coordinates [HMS05, Hor08]. The drawback of this method is that it 
depends on an a priori knowledge of separation coordinates and can therefore not 
be generalised. Besides, it is highly computational and thus unsuitable for higher 
dimensions. Nothing similar is known for the sphere or hyperbolic space. 

1.3. Method &: results. Our method is based on a purely algebraic description 
of the vector space of Killing tensors on constant curvature manifolds in combi- 
nation with the corresponding algebraic equivalent of the Nijenhuis integrability 
conditions. We consider the standard 3-sphere, isometrically embedded as the unit 
sphere in a 4-dimensional Euclidean vector space {V,g): 

^{v e V: g{v,v) = 1} c ^ . 

The vector space of Killing tensors K on is then naturally isomorphic to the 
vector space of algebraic curvature tensors on the ambient space V [MMS04] . These 
are (constant) tensors Raibiaob2 on V having the symmetries of a Riemannian cur- 
vature tensor. The isomorphism is explicitely given by 

K^{v,w) := Ra,b,a2b2x'''x'''v'''w''' xeS\ v,weT,S'^, (1.6) 

and equivariant with respect to the natural actions of the isometry group 0(V, g) on 
both spaces. Under this correspondence the Nijenhuis integrability conditions (1.4) 
for a Killing tensor are equivalent to the following algebraic integrability conditions 
for the associated algebraic curvature tensor [Schl2]: 

9r,R\a2b2R'd,C2d2=0 (l-7a) 
a.6.^V^^'..c... =0. (1.7b) 

Here the Young symmetrisers on the left hand side denote complete antisymmetrisa- 
tion in the indices 02 , 62 , C2 , ^2 respectively complete symmetrisation in the indices 
ai, 61, ci, di. These homogeneous algebraic equations define the space of integrable 
Killing tensors as a projective variety. 

The key to the solution of the algebraic integrability conditions is a reinterpre- 
tion in terms of the decomposition of algebraic curvature tensors into irreducible 
representations of the isometry group - the selfdual and anti-self dual Weyl part, 
the trace-free Ricci part and the scalar part. Under the action of the isometry 
group these components can be brought to a certain normal form, which simplifies 
considerably if the algebraic integrability conditions are imposed. 
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We regard an algebraic curvature tensor on y as a symmetric endomorphism on 
A'^V. In this sense we will say that an algebraic curvature tensor R is diagonal if 
Rijki = unless = {k,l}. In dimension four a diagonal algebraic curvature 

tensor can be parametrised as 

RoaOa — Wa + + Rp^Pf — Wa ^ + ]^ , (1-8) 

where (a,/?, 7) denotes any cyclic permutation of (1,2,3). In this parametrisation 
the Wa are the eigenvalues of the self-dual Weyl part, which turn out to be the 
same as those of the anti-self-dual Weyl part, the ta parametrise the eigenvalues of 
the trace free Ricci tensor and s is the scalar curvature. In a first step we prove 
that diagonal algebraic curvature tensors define a (local) slice for the action of the 
isometry group on integrable Killing tensors. 

Theorem 1.4 (Integrability implies diagonalisability). Under the action of the 
isometry group any integrable Killing tensor on S'^ is equivalent to one with a 
diagonal algebraic curvature tensor. 

The first algebraic integrability condition is identically satisfied on diagonal al- 
gebraic curvature tensors and the second reduces to a single algebraic equation on 
the diagonal entries. This yields a description of integrable Killing tensors with 
diagonal algebraic curvature tensor as a projective variety, equipped with the resi- 
dual action of the isometry group. We call this variety the Killing- Stackel variety 
(KS-variety). 

Theorem 1.5 (Integrability in the diagonal case). A Killing tensor on with 
diagonal algebraic curvature tensor (1.8) is integrable if and only if the associated 
matrix 

Ml -tg t2 \ 

h A2 -tl (1.9) 

\-t2 h A3 / 

has zero trace and zero determinant, where 

Ai W2 — A2 :— W3 — wi A3 :— Wi — W2 ■ 

This defines a projective variety in which carries a natural S^^-action, given by 
conjugating (1.9) with the symmetries of the regular octahedron in M? with vertices 
±ei and adjacent faces oppositely oriented. Two diagonal algebraic curvature ten- 
sors are equivalent under the isometry group if and only if their associated matrices 
(1.9) are equivalent under this S^-action. 

Combining Theorems 1.4 and 1.5 yields a solution to Problem II: The projec- 
tivisation of the quotient (1.5) is isomorphic to the quotient of the KS-variety by 
its natural S'4-action. Moreover, in the course of the proof of both theorems we 
establish isometry invariants for the integrability of a Killing tensor. 

Proposition 1.6 (Isometry invariants). Let W and T be the Weyl and Ricci part 
of the algebraic curvature tensor corresponding to a Killing tensor on C V, 
regarded as endomorphisms on A^V . Then this Killing tensor is integrable if and 
only if 

(i) tr[W,T]2 = 

(ii) tr (*W^) = tr (*W"^) = 0, where * denotes the Hodge operator on A^V and 
(Hi) the endomorphisms I, W and — T^ are linearly dependent. 
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With this result the check whether a given KiUing tensor on is integrable 
becomes algorithmically trivial. 

Naively, one can find the Stackel systems by solving the commutator equation 
[K, K] — for a fixed integrable Killing tensor K. But a far simpler way is to exploit 
the algebraic geometry of the KS-variety. Recall that Stackel systems are linear 
spaces of dimension n and always contain the metric. We show that the Killing 
tensors in a Stackel system have simultaneously diagonalisable algebraic curvature 
tensors. Therefore Stackel systems on correspond to certain projective lines in 
the KS-variety. In general, the projective lines in a projective variety constitute 
themselves a projective variety, called the Fano variety. That is, we are looking for 
a subvariety of the Fano variety of the KS-variety. By definition, the KS-variety 
is a linear section of the variety of 3 x 3 matrices with vanishing determinant, 
the so called full determinantal variety. Fortunately, the Fano variety of a full 
determinantal variety is well understood just in the case we need here. 

So we only have to use the commutator equation to check which projective lines 
correspond to Stackel systems, rather than to actually solve it. In order to do 
so, we derive an algebraic equivalent of the commutator equation in terms of the 
corresponding algebraic curvature tensors. This entails a detailed description of the 
algebraic geometry of the KS-variety and the Stackel systems therein. In particular, 
we can relate the algebraic geometric properties of the KS-variety to the geometric 
properties of the corresponding integrable Killing tensors. For this, we remark that 
there is a natural way how to extend (integrable) Killing tensors from S™ to 5'" 
for m < n. 

Theorem 1.7. For M — S'^ the projectivisation of the quotient (1.5) is isomorphic 
to the quotient of the KS-variety by its natural action of S/^. Furthermore: 

(i) The singularities of the KS-variety are extensions of the metric on S"^ re- 
spectively on . 

(ii) Stackel systems determine projective lines in the KS-variety which consist 

of matrices (1.9) annihilating a given vector. 
(Hi) Projective planes in the KS-variety that consist of matrices (1.9) annihi- 
lating a fixed vector correspond to extensions of (integrable) Killing tensors 
from 5*2 to . 

(iv) The projective plane of antisymmetric matrices in the KS-variety corre- 
sponds to the projective space of special Killing tensors. These are Killing 
tensors K of the form 

K = L-{tTL)g, 

where L is a special conformal Killing tensor^. 

Based on the algebraic description of Stackel systems we derive an algebraic de- 
scription of Benenti systems as well and answer Question I. The answer is "almost", 
in the following sense. On one hand, there are Stackel systems on S^, which are 
not Benenti systems. So the answer is clearly "no". But on the other hand these 
are exactly the extensions of Killing tensors from 5"^ and for S'^ the answer is easily 
seen to be "yes". 



cf. Definition 6.1 
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From the algebraic description of Benenti systems we derive the following para- 
metrisation of IC{S^), which is equivariant under the isometry group. This solves 
Problem I. 

Theorem 1.8 (Equivariant parametrisation). The algebraic variety JC{S^) is the 
closure of the set of algebraic curvature tensors of the form 

X2h©h + Xog®g ft. e Sym^(y), Ao,A2eM. 

Here ® denotes the Kulkarni-Nomizu product^ and ^g©g is the algebraic curvature 
tensor of the metric. That is, JC{S^) is the closure of the projective cone from the 
metric over the image of the Kulkarni-Nomizu square h i-^ h ® h, a map from 
symmetric tensors to algebraic curvature tensors. Remark that the apex lies itself 
in this image. 

Since the set of algebraic curvature tensors of the form h&h is invariant under the 
projective group PGL(F) of S^, the answer to Question II is "almost" as well. The 
reason is that although IC{S'^) is not projectively invariant, this "invariant part" 
still contains the essential part of each Stackel system unless it is an extension 
from S^. But again, on S"^ the answer is easily seen to be "yes". Our answers to 
Questions I and II indicate that one should think of the space of integrable Killing 
tensors on the sphere as composed of strata, where the non-generic strata come from 
lower dimensional spheres. Then each stratum is generated by Benenti systems and 
invariant under the projective group in the respective dimension. 

As a simple consequence of Theorem 1.7 we get a canonical representative in 
each Stackel system. 

Corollary 1.9. Each Stackel system on contains an essentially unique special 
Killing tensor. 

If the eigenvalues of the corresponding special conformal Killing tensor are sim- 
ple, they define the respective separation coordinates [CraOSa]. Moreover, from the 
so-called cone construction we know that a special conformal Killing tensor L is 
the restriction of a constant tensor L on the metric cone over the manifold [MM 10]. 
In our case this cone is just the ambient Euclidean space V. This is the reason 
why we can solve the eigenvalue problem for integrable Killing tensors explicitely. 
In this manner we get a purely algebraic classification of separation coordinates, 
which solves Problem IV. The details will be given later. We only mention that the 
natural parameters are the kernel of the matrix (1.9) and - in order to distinguish 
extensions - the eigenvalues of the trace free Ricci tensor. 

A solution to Problem III can be derived from the description of Stackel systems 
in the KS-variety. This gives the classification space for separation coordinates on 

modulo isometrics. 

Theorem 1.10 (Moduli space of separation coordinates). The space of Stackel 
systems modulo isometrics is homeomorphic to the quotient of the blow-up o/P^ in 
the four points (±1 : ±1 : ±1) € under the natural action of the symmetry group 
S4 of a tetrahedron inM.^ inscribed in a cube with vertices (±1,±1,±1). 

The fundamental domain of the S'4-action is a pentagon bounded by the non- 
generic separation coordinates. In fact, the multiplicities of the eigenvalues of the 



2cf. (6.4) 
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constant symmetric tensor L on V identify this pentagon to the associahedron K^, 
also known as StashefF polytope [Sta63, Lee89, Sta98]. This shows that the space of 
Stackel systems in the KS-variety is naturally homeomorphic to the permutoasso- 
ciahedron KP^ [Kap93]. These observations suggests a relation between separation 
coordinates on spheres and moduli spaces of stable curves of genus zero with marked 
points [DM69, Knu83]. In fact, the following can be proven [SV12J. 

Theorem 1.11. The space of Stackel systems on S"* with diagonal algebraic cur- 
vature tensors has a natural structure of an algebraic variety isomorphic to the real 
part A^o,n+2(]R) of the Deligne-Mumford-Knudsen compactification A4o_n+2{^) of 
the moduli space A^o.n+2(C) of stable algebraic curves of genus zero with n + 2 
marked points. 

1.4. Generalisations. Most of our results are stated for constant curvature man- 
ifolds of arbitrary dimension. What impedes a straightforward generalisation to 
higher dimensional spheres is that in the proof of diagonalisability we make essential 
use of the Hodge decomposition which only exists in dimension four. Theorem 1.11 
together with a proof of the following conjecture would show, that the classification 
space for separation coordinates on S" is the associahedron Kn+i- 

Conjecture 1.12. Under the action of the isometry group any integrable Killing 
tensor on S" is equivalent to one with a diagonal algebraic curvature tensor. 

In Lorentzian signature the Hodge star squares to minus one. A generalisation 
to hyperbolic space therefore naturally requires complexification. But over the 
complex numbers, all signatures are equivalent. This indicates that one should 
solve the algebraic integrability conditions over the complex numbers. The variety 
of complex solutions then contains not only the solution for hyperbolic space, but 
for all (non-flat) pseudo-Riemannian constant curvature spaces as real subvarieties, 
singled out by appropriate choices of a real structure. In dimension three there is 
an elegant way to do this using spinors and the Petrov classification of curvature 
tensors [MS 12]. The main complication comes from the fact that diagonalisability 
does not hold anymore in the complexified case. 

The biggest challenge remains without doubt a non-trivial classification of inte- 
grable Killing tensors and orthogonal separation coordinates on some non-constant 
curvature manifold. This is completely out of scope of common, purely geometric 
techniques. Promising candidates for the application of our algebraic approach are 
spaces of constant holomorphic sectional curvature, semi-simple Lie groups, homo- 
geneous spaces or Einstein manifolds. We propose a three-stage procedure, bearing 
in mind as a guiding example: 

(i) Identify the isomorphism class of the representation of the isometry group 
on Killing tensors and find an explicit isomorphism (as in [MMS04] for 
constant curvature manifolds). 

(ii) Translate the Nijenhuis integrability conditions for Killing tensors into 
purely algebraic integrability conditions (as in [Sclil2] for constant cur- 
vature manifolds). 

(iii) Solve the algebraic integrability conditions by means of representation the- 
ory and algebraic geometry (as in this article for S^). 

We consider the present work as a proof of concept that this algebraic approach to 
the problem of separation of variables is viable. 
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1.5. Structure of the article. The article is organised as follows: In Section 2 we 
examine some properties of algebraic curvature tensors that will be used to solve 
the algebraic integrability conditions in Section 3. We then describe the algebraic 
geometry of this solution in Section 4 and identify the Stackel systems therein in 
Section 5. In Section 6 we give various geometric constructions for integrable Killing 
tensors and interpret them algebraically. In the last section we demonstrate how 
the classification of separation coordinates can be recovered within our framework. 

Acknowledgements: I would like to thank Robert Milson for his motivation 
and numerous comments during inspiring discussions about my findings. He in 
particular spotted a missing equation that follows from the second algebraic inte- 
grability condition. I would also like to thank the organisers of the 36th University 
of Arkansas Spring Lecture Series "Conformal Differential Geometry and its Inter- 
action with Representation Theory" for their kind invitation to give a talk on my 
results. Finally I would like to thank Alexander P. Veselov for pointing out the link 
between my solution and moduli spaces of stable curves. 

2. Properties of algebraic curvature tensors 

Algebraic curvature tensors on a vector space V are four-fold covariant tensors 
Raibia2b2 having the same algebraic symmetries as a Rlemannlan curvature tensor: 



antisymmetry: ^&iai02h2 = —-^01610262 — ^01616202 (2.1a) 

pair symmetry: -^02620161 = ^01610262 (2-lb) 

Bianchi identity: Raxbta^b^ + Raia2b2bi + Raib2bia2 = (2.1c) 

Given a scalar product g on V, we can raise and lower indices. The symmetries 
(2.1a) and (2.1b) then allow us to regard an algebraic curvature tensor Raia2bib2 on 
F as a symmetric endomorphism R""^^^ a2bo '-'-'^ space K'^V of 2-forms on V . Since 
we will frequently change between both interpretations, we denote endomorphisms 
by the same letter in boldface. 

2.1. Decomposition. In the special case where dimV^ ~ 4, the Hodge star opera- 
tor " defines a decomposition 

of K^V into its ±1 eigenspaces. We can therefore write an algebraic curvature 
tensor R and the Hodge star as block matrices 

R = I I I + ^ I I " I * = ^^+ 1 " ^ (2.2a) 

= (2.2b) 
(2.2c) 
(2. 2d) 

Here J+ and /_ denote the identity on K^V respectively K'^V . The conditions 
(2.2b) assure symmetry, condition (2.2c) says that s = 2trR and condition (2. 2d) 
is a reformulation of the Blanchl Identity (2.1c). Indeed, the Bianchi identity is 



w+ 























+1+ 











with the 3 x 3-blocks satisfying 



tr VK+ + tr = 
tvW+-ivW- =0. 
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equivalent to the vanishing of the antisymmetrisation of i?aifcia2h2 ^our indices. 
In dimension four this can be written as s"'^''^"''^^^ Ra-^bia2b2 = or 

^a,2b2 -paibi n 

^ aibi-^ a2b2 — " ' 

where £'^2620161 jg ^j^g totally antisymmetric tensor. Notice that £°''^^^ciibi nothing 
else than the Hodge star operator and hence 

tr(*R)=0, (2.3) 

where tr : 'End{A'^V) ^- M is the usual trace. Now (2. 2d) is (2.3) applied to R in 
(2.2). 

The space of algebraic curvature tensors is an irreducible GL(y)-representation 
and (2.2) gives a decomposition of this representation into irreducible representa- 
tions of the subgroup SO{V,g) when dim!/ = 4. As the notation already suggests, 
we can relate these components to the familiar Ricci decomposition 

W + T + S 

of an algebraic curvature tensor R into 

• a scalar part 

*5'aifcia2fc2 ■ 12^.9 aia29b\b2 9aib29bia2^ 

given by the scalar curvature 

„ _ „aia2 6162 R . . 

• a trace free Ricci part 

Taib\a2b2 — 2 ("^"i '^2 3bl ^2 ~ T^aib29bia2 7fcia2.9aift2 "t' ^6162 3ai 02 ) 

given by the trace free Ricci tensor 

J^aia2 y J^aibia2b2 4i/aia2 

• and a Weyl part, given by the totally trace free Weyl tensor W R—T—S. 
It is not difficult to check that Tl =0 and hence 



W = 



w+ 
















T± 






s = 

















(2.4) 



We see that W+ and W- are the self-dual and anti-self-dual part of the Weyl tensor. 
Implicitly, the above interpretation also provides an isomorphism 

Sym„(y) A Hom(A^y,A2_y) 
T ^ T± 

between trace-free Ricci tensors and homomorphisms from self-dual to anti-self-dual 
2- forms. 
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2.2. The action of the isometry group. The Killing tensor equation (1.2) as 
well as the integrability conditions (1.4) are invariant under the action of the iso- 
metry group. In other words, the variety of integrable Killing tensors is invariant 
under isometrics. This will allow us later to put the algebraic curvature tensor 
of an integrable Killing tensor on to a certain normal form and to solve the 
algebraic integrability conditions for this particular form. Recall that the isometry 
group of S" C V is 0(V) and that SO(V^) is the subgroup of orientation preserving 
isometries. 

We first examine the induced action of SO(F) on algebraic curvature tensors of 
the form (2.2). The standard action of SO(F) on V induces a natural action on 
A'^V which is the adjoint action of SO(V^) on its Lie algebra under the isomorphism 

so{V) ^ A^V . 

Now consider the following commutative diagram of Lie group morphisms for 
dim F = 4: 



Spin(l^) 



-> Spin(A2_T/) X Spm{AlV) 



2:1 



2:1 



2:1 



tt: so(y) S0(A2_y) X SO{AlV) ^ SO{A^V) 

I [/_; 



u 



Here the double covering tt in the second row is induced from the exceptional 
isomorphism in the first row via the universal covering maps (vertical) . Under the 
induced isomorphism tt* of Lie algebras. 

so{V) so{AlV) e so{A^_V) , 

the adjoint action of an element U G SO(y) corresponds to the adjoint actions of 
U+ e SO(A^F) and C/_ e SO(Aiy). Since A\V has dimension three, the Hodge 
star operator gives an isomorphism 

so{AlV) ^ A^AlV ^ AlV 

under which the adjoint action of C/+ corresponds to the standard action on A'^V 
and similarly for [/_ . The Hodge decomposition completes the above isomorphisms 
to a commutative diagram 



so{V) 



A^V 



-> so{AlV) ® so{AlV) 



AlV 



AlV. 



which shows that the natural action of [/ e S0(1/) on A^V is given by 

fU: ' 



U- 



e 7r(S0(V")) C SO(AV) 



(2.6) 



Hence the natural action of U on End(A^y) is given by conjugation with this 
matrix. Restricting to algebraic curvature tensors we get: 
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Proposition 2.1 (Action of the isometry group). An element U E SO(V^) acts on 
algebraic curvature tensors in the form (2.2) by conjugation with (2.6), i. e. via 

W+ ^ UlW+U+ W- h-> U!^W^U- T± ^ U*_T±U+ , (2.7) 

where ■n{U) = {U^^U^) is the image of U under the double covering 

tt: S0{V) ^ SO{h.\V) X SO(A^y) . 

We can describe the twofold cover it expHcitely in terms of orthonormal bases on 
the spaces V , K^jV and h^JV . It maps an orthonormal basis (cq, ei, 62, 63) of V to 
the orthonormal bases (ry+i, 77+2, '7+3) of A^T^ and (ry_i, ?7_2, ry-s) of A^F, defined 
by 

jyia := ^(eo A ea ± 6/3 A e-y) (2.8) 

for each cyclic permutation (q;,/3,7) of (1,2,3). From this description also follows 
that kerTT = {±/}. 

The action of 0{V) on algebraic curvature tensors is now determined by the 
action of some orientation reversing element in OiV), say the one given by reversing 
the sign of eg and preserving ei, 62 and 63. This element maps •q±ci to —rjz^a- Hence 
its action on algebraic curvature tensors in the form (2.2) is given by conjugation 
with 

respectively by mapping 

W+ T±^T^. (2.9) 

2.3. Aligned algebraic curvature tensors. W+ and W- are symmetric and 
hence simultaneously diagonalisable under the action (2.7) of SO(T^). On the other 
hand, the singular value decomposition shows that T± is also diagonalisable under 
this action, although in general not simultaneously with W+ and . The following 
lemma gives a criterion when this is the case. 

Lemma 2.2. Let W+ and W- be symmetric endomorphisms on two arbitrary Eu- 
clidean vector spaces A_|_ respectively A_ and suppose the linear map T± : A_|_ — )■ A_ 
satisfies 

T±W+ = W-T± . 

Then there exist orthonormal bases for A+ and A_ such that W+, W- and T± are 
simultaneously diagonal with respect to these bases. 

Proof. It suffices to show that we can chose diagonal bases for W+ and W- such 
that the matrix of T± has at most one non-zero element in each row and in each 
column, for the desired result can then be obtained by an appropriate permutation 
of the basis elements. The above condition implies that T± maps eigenspaces of 
W+ to eigenspaces of W- with the same eigenvalue. Without loss of generality we 
can thus assume that A+ and A_ are eigenspaces of W+ respectively VF_ with the 
same eigenvalue. But then W+ and W- are each proportional to the identity and 
therefore diagonal in any basis. In this case the lemma follows from the singular 
value decomposition of T±. □ 

Lemma 2.3. The following conditions are equivalent for an algebraic curvature 
tensor (2.2).' 
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(i) T±W+ = W-T± 
(li) [W, T] = 

(iii) W and T are simultaneously diagonalisable under SO{A^V). 

Proof. The equivalence of (i) and (ii) follows from (2.4). For the equivalence of 
{a) and (iii) it suffices to note that W and T are both symmetric and hence 
diagonalisable. □ 

The following terminology is borrowed from general relativity. 

Definition 2.4. We say that an algebraic curvature tensor is aligned, if it satisfies 
one of the equivalent conditions in Lemma 2.3. 

By Lemma 2.3, for any aligned algebraic curvature tensor we can find an or- 
thonormal basis of V such that 



\ u;+2 I \ u;_2 | T± = | ^2 ] (2.10a) 

with 

w+i + u)+2 + w+3 = W-i + u)_2 + W-3 = . (2.10b) 

To simplify notation we agree that henceforth the indices (a, /3,7) will stand for 
an arbitrary cyclic permutation of (1,2,3). Changing the basis in A^V from (2.8) 
to A Cj, ^ i < J ^ 4, we obtain the independent components of an aligned 
algebraic curvature tensor: 

7/J ( „ 4- 7/1 „ 

-\-t + — 







^ ti \ 


1-1 


t2 




^ <3 / 







-« + 






2 




RjS-iP-i = 




-a + 






2 








-a 






2 





to. + ^ (2.11) 



The trace free Ricci tensor of this aligned algebraic curvature tensor is diagonal 
and given by 

Too = ta + tp + tj ^ Too + Taa ,„ , „\ 

T "~ 2 ■ ^ 

This is nothing else than the restriction of the isomorphism (2.5) to diagonal tensors. 
2.4. Diagonal algebraic curvature tensors. 

Definition 2.5. We call an algebraic curvature tensor on V diagonal in a basis 
Ci if it is diagonal as an element of End(A^T^) with respect to the associated basis 
Ci /\ Cj, i < j. We call it diagonalisable if it is diagonalisable under the adjoint 
action of SO{V) on A^V, i.e. under the subgroup n{SO{V)) C SO{A'^V). 

Of course, being a symmetric endomorphism on A'^V, an algebraic curvature 
tensor is always diagonalisable under the full group SO(A^V^). 

Proposition 2.6 (Diagonalisability criterion). An algebraic curvature tensor is 
diagonalisable if and only if it is aligned and W+ has the same characteristic poly- 
nomial as W- . 
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Proof. Consider an aligned algebraic curvature tensor and suppose that W+ has the 
same characteristic polynomial as Then we can find a transformation (2.7) 

such that, as 3 X 3 matrices, W+ = W-. The condition W-Tz^z — T±W+ then 
implies that T± can be simultaneously diagonalised with W-^- = W-. This means 
we can assume without loss of generality that 

W+a = W-a Wa 

for a = 1,2,3 in (2.10), so that (2.11) reads 

RoaOa ~ Wa + ta + Roap-y — RfSjPj — Wa — + (2.13a) 

with Wa subject to 

Wi+W2 + w^ = Q. (2.13b) 

Obviously this algebraic curvature tensor is diagonal. On the other hand, any 
diagonal algebraic curvature tensor can be written in the form (2.13), because 

s — 2 (i?oioi + -R2323 + -R0202 + -R3131 + ^0303 + -R1212) 

, _ RoaOa — Rp^lSj 
ta - ^ 

-RoqOq + Rp-f/}-/ s 

= 2 12- 

□ 



(2.14a) 
(2.14b) 

(2.14c) 



2.5. The residual action of the isometry. For later use we need the stabiliser 
of the space of diagonal algebraic curvature tensors under the action of the isometry 
group. Such tensors are of the form (2.10) with w^a = and are invariant under 
the transformation (2.9). Hence it is sufficient to consider only the subgroup SO(y) 
of orientation preserving isometrics. 

Let U e SO(F) with Tr{U) = (U+,U^) be an element preserving the space of 
diagonal algebraic curvature tensors under the action (2.7). Then U+ G SO{A^V) 
and U- G SO{A^V) preserve the space of diagonal matrices on A'^V respectively 
A^y. In the orthogonal group, the stabiliser of the space of diagonal matrices 
under conjugation is the subgroup of signed permutation matrices. They act by 
permuting the diagonal elements in disregard of the signs. In particular, in S0(3) 
the stabiliser subgroup of the space of diagonal matrices on is isomorphic to ^4. 
Under the isomorphism 6*4 = ^3 K V4, the permutation group S'3 is the subgroup 
of (unsigned) permutation matrices and acts by permuting the diagonal elements, 
whereas the Klein four group V4 is the subgroup of diagonal matrices in SO (3) and 
acts trivially. 

Using the above, it is not difficult to see that the stabiliser of the space of diagonal 
algebraic curvature tensors under the action (2.7) is isomorphic to 5*4 x kerTr, with 
kerTT = {±/} acting trivially. Under the isomorphism 6*4 = S'3 x V4, the factor S'3 
acts by simultaneous permutations of (wi, 102,103) and (^1,^21^3) and the factor V4 
by simultaneous flips of two signs in (^1,^2,^3)- 
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3. Solution of the algebraic integrability conditions 
This section is dedicated to proving Theorems 1.4 and 1.5. 



3.1. Reformulation of the first integrability condition. In the same way as 
for the Bianchi identity we can reformulate the first integrabihty condition (1.7a) 
as 

>C2d2 



£12 02 C2a2 







where 



r(R*R) =0, 
End(A2F) — > End(F) 



(3.1) 



(3.2) 



kl ^ kj 

denotes the Ricci contraction. Notice that E := R * R is symmetric, but does in 
general not satisfy the Bianclii identity. A proof of the following lemma can be 
found in [ST69]. 

Lemma 3.1. The kernel of the Ricci contraction (3.2) is composed of endomor- 
phisms E satisfying 

*E* = E* tr E = . 

Applying this to the symmetric endomorphism 
/ 

E := R* R 

V 



{W+ + ^I+r - T^T± 




_T^(W^_ + ^/_) 


T±{W+ + f-^I+) 

-{W^ + ^/_)T± 


T±T^ 


- (T^_ + ^/_)2 



in (3.1) and using (2.2c) we get 



tr Wi = tr 



(3.3a) 
(3.3b) 



This shows: 



Proposition 3.2. The first algebraic integrability condition for an algebraic cur- 
vature tensor R ~ W + T with Weyl part W and (not necessarily trace free) Ricci 
part T is equivalent to 



[W,T] = 



tr(*W^ 



0. 



In particular, an algebraic curvature tensor satisfying the first algebraic integrability 
condition is aligned. 

As a consequence of Proposition 2.6 we get: 

Corollary 3.3. A diagonalisable algebraic curvature tensor satisfies the first alge- 
braic integrability condition. 
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3.2. Integrability implies diagonalisability. The aim of this subsection is to 
prove Theorem 1.4. In view of Propositions 2.6 and 3.2 it is sufficient to show that 
if an algebraic curvature tensor satisfies the algebraic integrability conditions, then 
W+ and W- have the same characteristic polynomial. We will first prove that the 
first algebraic integrability condition implies that W+ has the same characteristic 
polynomial either as +W- or as —W-. We then prove that the latter contradicts 
the second algebraic integrability condition. 

We saw that an algebraic curvature tensor which satisfies the first algebraic 
integrability condition is aligned and thus of the form (2.10) in a suitable orthogonal 
basis. The first integrability condition in the form (3.3) then translates to 

w'^^ + + ■"'+3 = ^'^ '"'-I + ^-2 + ^-3 — ^'^ (3.4a) 

for some r ^ and 

w+ata = taW^a ■ (3.4b) 

If we regard (w+i, w-|_2, w+3) and (w-i, i(;_2, 1^-3) as vectors in M'^, then each equa- 
tion in (2.10b) describes the plane through the origin with normal (1,1,1) and each 
equation in (3.4a) the sphere of radius r centered at the origin. Hence the solutions 
to both equations lie on a circle and can be parametrised in polar coordinates by a 
radius r ^ and angles ip+ , if- as 

w+a ~ r cos[(p+ + a^) w_q = r cos((/3_ + a^) . (3.5) 

Assume now that W+ and W- have different eigenvalues. Then we have T± = 0. 
Indeed, if 7^ for some a, then (3.4b) shows that io+q = W-a for this index 
a. But then (3.5) implies that W+ and have the same eigenvalues, which 
contradicts our assumption. Hence T± = 0. 

Without loss of generality we will set s = in what follows, because integrabi- 
lity is independent of the scalar curvature s. This follows from the fact that the 
Killing tensor corresponding to the scalar component in the Ricci decomposition is 
proportional to the metric and thus to the identity when regarded as an endomor- 
phism. But adding a multiple of the identity only alters the eigenvalues, not the 
eigendirections, and therefore has no effect on integrability. 

After these considerations, the independent components (2.11) of our algebraic 
curvature tensor are 

U - + - U O _ ■^+" ^ R\ 
noaOa — ^ — -tipflS-) noafj^ — ^ . W-DJ 

Choosing ai = 61 = ci = di = in the second integrability condition (1.7b) for 
this algebraic curvature tensor yields 



— RiQa2b2^jOkoRlOc2d2 — 

C2 ' ^ 

"2 



~ -RaOa2b2^aOaoRaOc2d2 — , 



since only terms with i— j = k = l^O contribute to the contraction. For 
(02,62,02,(^2) = (0,1,2,3) this is gives 

RaOOaRaOaoRaOI3'-f = {RoaOa ) ^ Rpa/S-y — 0, 

where the sums run over the three cyclic permutations {a,f3,^) of (1,2,3). We 
conclude that 

^ {w+a + W^af{w+a - W-a) ^ ■ 
a=l,2,3 
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Substituting (3.5) into this equation yields, after some trigonometry, 



3 „„2 V'+ - "P- „:„ ¥'+ - f - „;„ 3 , 



r cos ^ sm ^ sm -(^+ + = . 

Substituting the solutions 

r = = + fcTT (p+ = — (/?- + /c^ G Z 

of this equation back into (3.5) shows that the set {w-i, W-27 ^-3} is either equal to 
{w+i, u'+2, or to {— w+i, — ^w+s}- The first case is excluded by assump- 
tion and we conclude that W+ and — PF_ have the same eigenvalues. Consequently 
we can choose an orthonormal basis of V such that 

W+Q = -W-a =: Wa ■ (3.7) 

Then the only remaining independent components of the algebraic curvature tensor 
(3.6) are 

-Roq/37 = Wa ■ (3.8) 

Now consider the second integrability condition (1.7b) for (02, 62, C2, d-i) — (0, 1, 2, 3), 
written as 

VA^yVMiW ^ 9^^ ^ih^ld^^'ia^kc^ = (3.9) 

with 

1 H 

Eibildi ■= -^-^ Ribia2b2RldiC2d2 (3.10) 

The product Ribia2b2Ridxc2d2 is zero unless {i, 61,02, 62} = C2, ^2} = {0, 1, 2, 3}. 

Since {02, 62, C2, ^2} = {0, 1, 2, 3}, we have Eib^id^ = unless {z, bi,l, di} = {0, 1, 2, 3}. 
The symmetries 

Ebiiidi — —Eibiidi = Eib^dii Eibj^idi = Euhibi (3-11) 

then imply that the only independent components of (3.10) are 

= < ■ (3-12) 
Now consider the full contraction of (3.9): 

Since the Ricci tensor of (3.8) is zero, this is results in 

E'''''R,jki + E'^'''Rak,^Q- 

Using the symmetries (3.11) of Eijki and those of Rijki, we can always permute a 
particular index to the first position: 

And since both iJ'-''^' and Rijki vanish unless k, 1} ~ {0, 1, 2, 3}, we get 

y^(-5'"°^^-Roa^7 + E^'^'^°'Ro^pa) + ''^{E^"'^'^ Rof Pa + E^^^"" R^ap-^) = . 

Here again, the sums run over the three cyclic permutations (a,/?, 7) of (1,2,3). 
Substituting (3.8) and (3.12) yields 

^{wl + wi^) - ^{wiw^ + wl^Wa) = . 

With (2.13b) the second term can be transformed to 

^(w^W^ + W^Wq) = ^(w^Wq -I- w'pW^) = ^ wfjiWa + W7) = - ^ , 
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implying 



By (2.13b) and Newton's indentity 



wl + W2+ wl 



3wiW2W3 — 3{wiW2 + W2W3 + W3Wi){wi + W2 + W3) + {wi + W2 + ICa)^ 



this is equivalent to 



W1W2W3 — . 



Using (2.13b) once more, this shows that {wi,W2,W3} = {—w,Q,+w} for some 
w S K. But then (3.7) implies that {w+i, u;+2, w+a} — {w-i,W-2tW-3}, which 
contradicts our assumption that the eigenvalues of W+ and W_ are different. Con- 
sequently W+ and W- have the same eigenvalues. This finishes the proof of Theo- 
rem 1.4: Any integrable Killing tensor on has a diagonalisable algebraic curva- 
ture tensor. 

3.3. Solution of the second integrability condition. It remains to prove Theo- 
rem 1.5. In the preceding subsections we showed that an algebraic curvature tensor 
which satisfies the algebraic integrability conditions is diagonal in a suitable ortho- 
gonal basis of V and that any diagonal algebraic curvature tensor satisfies the first 
of the two algebraic integrability conditions. We therefore now solve the second 
integrability condition (1.7b) for a diagonal algebraic curvature tensor in the form 
(2.13). To this aim consider the tensor 



appearing on the left hand side of (1.7b). Suppose it does not vanish. Then 
we have bi G {02,62} and di G {02,^2}, so without loss of generality we can 
assume bi = 62 and di = c?2- Then only terms with i = j = a2 and k = I = C2 
contribute to the contraction and hence {01,02} — {ci,C2}. If now 02 = 02, then 
(3.13) vanishes under complete antisymmetrisation in 02, 62: C2, ^2- This means that 
the left hand side of (1.7b) vanishes unless {oi, &i, ci, di} — {02,62,22,^2}- Note 
that 02, 62, C2, d2 and therefore oi, 61, ci, di have to be pairwise different due to the 
antisymmetrisation. 

For dimy = 4 in particular, (1.7b) reduces to a sole condition, which can be 
written as 



(3.13) 



/ 1 




det 



1 



V 1 

Substituting (2.13) yields 



det 




(3.14) 



The form of this equation allows us to give an isometry invariant reformulation of the 
algebraic integrability conditions, which does not rely on diagonalising the algebraic 
curvature tensor. This also entails the isometry invariants in Propositions 1.6. 
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Proposition 3.4. The algebraic integrahility conditions for an algebraic curvature 
tensor with WeylparfW and (not necessarily trace free) RiccipartT are equivalent 
to 

[W,T]=0 (3.15a) 

tr(*W2) ^ tr(*W3) (3.15b) 

and the linear dependence o/ W, — and the identity matrix. 

Proof. This is a consequence of Theorem 1.4 and Proposition 2.6, together with 
Corollary 3.3 and the second integrahility condition in the form (3.14). Indeed, 
(3.15a) is the definition of alignedness for R and tr(*W''') = is equivalent to 
trW^ = tiWll. For fc = 1 this equation is the Bianchi identity in the form (2. 2d). 
Hence (3.15b) means that W+ and W- have the same characteristic polynomial. 
The linear dependence of W, — and the identity matrix is equivalent to 
(3.14). □ 

For ti — t2 — t^ ^ the left hand side of (3.14) is the Vandermonde determinant. 
Therefore 

{Wi - W2) tl + {W2 - W3) t\ + {W3 -wi)tl 

+ (Wi - W2)(W2 - W3)(W3 - Wl) = . (3.16) 

This proves the first part of Theorem 1.5, for this equation can be written in the 
form 

I'w2-W3 ~t3 <2 \ 

0. 

The second part follows from the residual action of the isometry on diagonal alge- 
braic curvature tensors as described in Subsection 2.5. 

Regarding the eigenvalues wi, W2, W3 S K as parameters subject to the restriction 

W1+W2 + W3 — O, (3.17) 

we can solve this equation for ii,i2,^3 G depending on the number of equal 
eigenvalues: 

(i) Ifwi, W2, W3 are pairwise different with < wp < w-y then (3.16) describes 
a one sheeted hyperboloid 

tiy \ ( tp \ I 




an 



1 



with semi axes 



«a = A 






- Wa) 


ap = 1 






-Wp) 


= 1 









ill) If Wa — Wp ^ Wry. then (3.17) implies that the solutions of (3.16) are those 
(^1,^2,^3) € with tp = zLta. Since (ta,—tp,—t^) and (ta^tp.t.^) are 
equivalent under the isometry group, we can assume tp = t^ without loss 
of generality. 
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(iii) If wi = W2 — W3, then (3.17) implies that wi = W2 = W3 = 0. In this 
case (3.16) is satisfied for any choice of ti,t2,t^ G M. In other words, an 
algebraic curvature tensor with zero Weyl component is always integrable. 

4. The algebraic geometry of the KS-variety 

We have seen that under isometrics any integrable Killing tensor is equivalent to 
one whose algebraic curvature tensor is diagonal and that the space of such tensors 
is described by the following variety. 

Definition 4.1. Let A = MJ' be the vector space of matrices whose symmetric part 
is diagonal and trace free, i. e. of matrices of the form 







-h 






M = 




A2 




(4.1) 








Ag^ 





with 

tr M = . 

We denote hy V <Z A the algebraic variety defined by the equation 

det M = 

and call its projectivisation P V C P.4 ^ P*^ the Killing- Stdckel variety (KS-variety). 
A matrix in V will be called a Killing-Stdckel matrix (KS-matrix). 

Remark 4.2. The residual isometry group action on diagonal algebraic curvature 
tensors defines a natural action of the permutation group Si on the Killing-Stackel 
variety. This action is given by conjugation with matrices in SO (3) under the 
embedding 6*4 C S0(3) defined in Subsection 2.5. Later it will be useful to consider 
Si C SO (3) as the symmetry group of an octahedron in M'^ with vertices ±6^ and 
adjacent faces opposedly oriented. 

The Killing-Stackel variety P V is a 3-dimensional projective subvariety in P^ and 
the quotient P V /S'4 encodes all information on integrable Killing tensors modulo 
isometries, the metric and scalar multiples. In this section we will investigate its 
structure from a purely algebraic geometric point of view. This will finally lead to 
a complete and explicit algebraic description of Stackel systems. 

Recall that Stackel systems are n-dimensional vector spaces of mutually commu- 
ting integrable Killing tensors. Note that every Stackel system contains the metric 
and that the metric corresponds to the zero KS-matrix. Therefore, if we assume 
that the algebraic curvature tensors of the Killing tensors in a Stackel system are 
mutually diagonalisable, then Stackel systems on S^ correspond to 2-dimensional 
planes in V or projective lines in the projectivisation P V. We will see a posteriori 
that our assumption is justified. Hence Stackel systems constitute a subvariety of 
the variety of projective lines on P V, also called the Fano variety of P V and denoted 
by Fi(PV). 

Let 2?3 be the generic determinant al variety. This is the variety of 3 x 3 matrices 
with vanishing determinant or, equivalently, matrices of rank one or two. By defi- 
nition P V is a subvariety of P2?3. Hence its Fano variety Fi(P V) is a subvariety of 
the Fano variety FiiWD^). But the latter is well understood. It contains different 
types of projective hnear spaces: 
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(i) projective spaces of matrices whose kernel contains a fixed line, 

(ii) projective spaces of matrices whose image is contained in a fixed plane and 

(iii) the projective plane of antisymmetric matrices. 

This motivates why we seek for such spaces in PV. Note that the projective plane 
of antisymmetric matrices is obviously contained in PV. 

Definition 4.3. We call a projective subspace in the KS-variety consisting of ma- 
trices whose kernel contains a fixed line (respectively whose image is contained in 
a fixed plane) an isokernel space (respectively an isoimage space). 

A 3 X 3 matrix M of rank two has a 2-dimensional image and a 1-dimensional 
kernel. Both are given by its adjugate matrix Adj M, which is the transpose of the 
cofactor matrix and satisfies 



(Adj M)M = Af(Adj M) = (detM)/ 
The adjugate of a matrix M e V is 

Adj M = tit2 ~ Asts t2t2 
\ tit3 + A2t2 t2h 



Ash 

A3A1 

Ai^i 



ts^i-Asis 

t3t2 + Aiti 

tata + AiAa 



(4.2) 



(4.3) 



The adjugate matrix of a rank two 3x3 matrix has rank one, which means that 
the columns of Adj M span a line in M.^. And because det Af = 0, we deduce from 
(4.2) that the kernel of M is the column space of Adj M. 

The adjugate matrix of a rank one 3x3 matrix is zero. It is known, that 
the singular locus of P2?3 is the subvariety of rank one matrices. The following 
proposition characterises the singular locus of the KS-variety. For simplicity we 
will be a bit imprecise and not strictly distinguish between a matrix M S V \{0} 
and its image KM G PV under the canonical projection V\{0} — > PV. Depending 
on the context, we will refer to an element of P V as "matrix" or "point". 



Proposition 4.4 (Singularities of the KS-variety). The KS-variety 
singular points, given by the six rank one matrices 



'V has ten 





/O 






f-l 


±l\ 




f+1 






+1 




V±2 ■■ = 








V±3 


±1 


-1 




\0 ±1 






Wi 


+1/ 




\o 






and the four skew symmetric matrices 
Co := V+i + V+2 + V+3 



V+ 



+7 ' 



(4.4a) 



(4.4b) 



where (a,/?, 7) denotes an even permutation of (1,2,3). Moreover, the rank one 
matrices m PV are exactly those six in (4.4a). 

Definition 4.5. We will call the six singularities V±a the rank one singular points 
in PV and the four singularities C; the skew symmetric singular points in P V. 

Proof. P V is singular at MA/ if and only if V is singular at Af . Since V is the zero 
locus of the determinant function det : ^ — > E, this is the case if and only if the 
derivative of the determinant function at Af G V, given by 



(ddet) 



M ■ 



A 

A 



tr(A Adj M) , 



(4.5) 
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is the zero map. The condition that tr(v4 Adj M) ~ for all A e ^ is equivalent to 

tl + A2A3 ^tl+ A3A1 ^tl + Ai A2 Alii = A2i2 Agtg = . 

We leave it to the reader to verify that the solutions M e V of these equations are 
exactly the matrices (4.4) and their multiples. 

The last statement now follows from (4.5) and the equivalence of rankAf — 1 
and Adj M = for 3 x 3-matrices. □ 

Remark 4.6. Under the natural S'4-action on PV, all rank one singular points are 
equivalent. This means that the singularieties V±a are all mapped to a single point 
V in the quotient PV /S'4. The same holds for the skew symmetric singular points, 
which are mapped to a single point C G P V /S'4 . 

We now compute the isokernel spaces in PV. To find all matrices Af G V 
annihilating a given vector n e M'^ , we consider the equation Mn = as an equation 
in M for a fixed ft = [ni, n2,n^), where M is of the form A/ = A + T with 

Ml \ 
A= A2 T = 

yo Aa/ 

We can regard t = (^1,^2 7^3) € K'^ as a parameter and solve An = —Tn for A. 
For the time being, let us assume nin2n3 ^ 0. We then obtain a linear family of 
matrices 

/tt3-t^2 -^3 h 

h ^ii - ^t3 1 . (4.6a) 




M 



\ -t2 h ^t2 ~ 



"3 ^ "3 

All these matrices satisfy det M = and the condition tr A^f = imposes the 
restriction 

^.^/n3_n2 ni_n3 n2_M 
V"2 713 ns Til ni 712/ 

on the parameter i S M^. This shows that the matrices in V annihilating a fixed 
vector n 7^ form a linear subspace. If the right hand side of (4.6b) is zero, i. e. 
if the condition is void, this subspace has dimension three and defines a projective 
plane in PV. This happens if and only if 

Iml = |«2| - l^al- (4.7) 

If the right hand side of (4.6b) is not zero, this subspace has dimension two and 
defines a projective line in P V. One can check that this also holds true for the case 
711712713 = 0. We have shown the following: 

Lemma 4.7. If n G M.^ satisfies (4.7), then the set of matrices m PV annihilating 
n is a projective plane of the form (4.6a). Otherwise it is a projective line. 

Recall that the kernel of a rank two matrix of the form (4.1) is the row space of 
the adjugate matrix (4.3). This defines a rational map 

tt: PV — ^ P2 



whose fibres are the maximal isokernel spaces. This map is well defined except for 
the six rank one matrices in PV. 
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We want to give a parametrisation of the isokernel lines that is uniform and more 
geometric than (4.6). For this purpose we define two embeddings l^v: PV 
given by 



for n — (ni : n2 : n^) E P^. One readily checks that L{n) and i'{n) indeed lie in 
P V. Note that the image of t is the projective plane of skew symmetric matrices 
in P V. Under the inclusion P V C PA = P^ the map t is a linear embedding of P^ 
in P^ whilst ly is the composition of the Veronese embedding 



and a projection P^ P^ under which the Veronese surface remains smooth. 

Proposition 4.8 (Isokernel spaces in the KS-variety). Recall from Proposition 4-4 
that the singular points in the KS-variety P V are the four skew symmetric singular 
points Co,Ci,C2,C3 as well as six rank one singular points V±a, a = 1,2,3. 

(i) The projective plane through V+i, V+2 o,nd V+3 is an isokernel plane in PV 
and containes Cq. In the same way, the three points V-a, V+ii o,nd 
define an isokernel plane through Cq. for a — 1^2, 3. 

Let the maps tt: P V P^ and t, i^: P'^ — > P V he defined as above. 

(ii) For any non-singular point M G P V, the points M , L{ir{M)) and z^(7r(Af)) 
are collinear hut do not all three coincide. The projective line they define lies 
m PV and is an isokernel space through M. Each of these lines cointains 
a unique skew symmetric point, namely l{-k{M)). 

Moreover, any isokernel space m P V of maximal dimension is of either of the above 
forms and hence contains a unique skew symmetric point. 

Remark 4.9. A corresponding characterisation of isoimage spaces follows from the 
fact that matrix transposition defines an involution of P V which interchanges iso- 
kernel and isoimage spaces. 

Proof, (i) follows directly from Lemma 4.7 and the explicit form (4.4) of the singular 
points. For (ii) notice that n :— 7r(M) is the kernel of M and well defined for 
M ^Va. The definitions (4.8) of t and v then show that n is also the kernel of t(n) 
and v{n). These definitions also show that M, i{n) and v{n) do not all coincide 
for M 7^ Ci. Their coUinearity follows from Lemma 4.7 in the case where n does 
not satisfy (4.7) and from t(ri) — v{n) in case it does. The last statement is now a 
consequence of Lemma 4.7. □ 

We can illustrate the content of Proposition 4.8 geometrically as follows. This 
is depicted in Figure 1 . The six rank one singular points of P V C P* constitute 
the vertices of a regular octahedron in P** whose faces (and the planes they define) 
are entirely contained in P V. The set of eight faces is divided into two sets of non- 
adjacent faces, corresponding to the four isokernel planes (shaded) respectively the 
four isoimage planes in PV. Opposite faces intersect in their respective centres. 
These are the skew symmetric singular points of P V and the intersection points of 




(4.8) 



p2 ^ p5 

(ni : ^2 : ^3) ^ [n\ : n\ : n1 : n2n^ : ^3711 : 7^1712) 
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the octahedron with the projective plane of skew symmetric matrices. In Figure 1 
this ninth projective plane in P V is depicted as the insphere of the octahedron. 

^+1 




Figure 1. Singularities and projective planes in the KS- variety. 

The permutation group S4 acts on this configuration by symmetries of the oc- 
tahedron preserving the two sets of non-adjacent faces. Matrix transposition cor- 
responds to a point reflection, which exchanges opposite faces and completes the 
S'4-action to the full octahedral symmetry. 

In anticipation of the results in the next section we introduce the following name: 

Definition 4.10. We call the isokernel lines of the form (ii) "Stdckel lines". 

Every Stackel line intersects the projective plane of skew symmetric matrices 
transversely. Conversely, every point on this plane determines a unique Stackel 
line unless it is one of the four singular skew symmetric points. The Stackel lines 
through these points are exactly those projective lines in each isokernel face which 
pass through the respective face center. 

5. Stackel systems 

In this section we express the condition that two Killing tensors commute as 
an algebraic condition on the corresponding algebraic curvature tensors. This will 
eventually justify the term "Stackel lines" chosen in Definition 4.10 by proving that 
these lines indeed correspond to Stackel systems. 

Remark 5.1. Most results in this and the following section are valid for an arbitrary 
non-fiat constant curvature (pseudo-)Riemannian manifold M that is embedded 
into some vector space V as the hypersurface 

M = {v eV: g{v,v) = 1} C V 
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of length one vectors with respect to a non-degenerated symmetric biHnear form g 
on V . Note that any non-flat constant curvature (psudo-)Riemannian manifold is 
locally isometric to such a model. 

5.1. Young tableaux. Throughout this section we use Young tableaux in order 
to characterise index symmetries of tensors. Young tableaux define elements in the 
group algebra of the permutation group Sd- That is, a Young tableau stands for a 
(formal) linear combination of permutations of d objects. In our case, these objects 
will be certain tensor indices. For the sake of simplicity of notation we will identify 
a Young tableau with the group algebra element it defines. It will be sufficient for 
our purposes to consider the examples below. For more examples and a detailed 
explanation of the notation as well as the techniques used here, we refer the reader 
to [Schl2]. 

A Young tableau consisting of a single row denotes the sum of all permutations 
of the indices in this row. For example, using cycle notation, 

\a2W\c2\ ^ e + (a2Ci) + (C1C2) + (c2a2) + (a2CiC2) + (020102) . 

This is an element in the group algebra of the group of permutations of the indices 
02, ci and C2 (or any superset). In the same way a Young tableau consisting of 
a single column denotes the signed sum of all permutations of the indices in this 
column, the sign being the sign of the permutation. For example. 



(ai5i) - (61^2) - (^201) + (0161^2) + (^26101) 



We call these row symmetrisers respectively column antisymmetrisers. The reason 
we define them without the usual normalisation factors is that then all numeri- 
cal constants appear explicitly in our computations (although irrelevant for our 
concerns) . 

The group multiplication extends linearly to a natural product in the group 
algebra. A general Young tableau is then simply the product of all column an- 
tisymmetrisers and all row symmetrisers of the tableau. We will only deal with 
Young tableaux having a "hook shape", such as the following: 



01 


02 Cl C2 1 


bi 




d2_ 





\al\a2\cl\c2\ 



(5.1) 



The inversion of group elements extends linearly to an involution of the group 
algebra. If we consider elements in the group algebra as linear operators on the 
group algebra itself, this involution is the adjoint with respect to the natural in- 
ner product on the group algebra, given by defining the group elements to be an 
orthonormal basis. Since this operation leaves symmetrisers and antisymmetrisers 
invariant, it simply exchanges the order of symmetrisers and antisymmetrisers in a 
Young tableau. The adjoint of (5.1) for example is 





02 1 Cl C2 


61 









— 

d2 



(5.2) 
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Properly scaled, Young tableaux define projectors onto irreducible Sd-representations. 
We are now able to state a lemma which expresses the decomposition 



given by the Littlewood- Richardson rule, in terms of Young projectors. Its proof 
can be found in [Schl2]. 

Lemma 5.2. In the group algebra of the permutation group of the indices ai, 02, 
Ci, C2, 61 and ^2 the following identity holds: 



1 
2! 



1 

4! 



10368 



ai 


02 Cl C2 1 




12 Cl C2 


h 








d2_ 




d2_ 





34560 



ai|n2|Ci|C2 



"1 |a2 I Ci I C2 I 



(5.3) 



The action of the permutation group on tensors extends linearly to an action of 
the group algebra. In particular, any Young tableau acts on tensors with corre- 
sponding indices. For example. 



T, 



a2h\C2 



TbiC2a2 ^^£10262 TaiC2b2 + 2^cih2a2 



To give another example, the operator (5.1) acts on a tensor Ti)-^i,2did2a2C2 by an 
antisymmetrisation in the indices 61 , 02 , C2 and a subsequent symmetrisation in the 
indices 61, 621 di, d2- In the same way its adjoint (5.2) acts by first symmetrising 
and then antisymmetrising. 

5.2. Commuting Killing tensors. In the following we express the fact that two 
Killing tensors commute (as endomorphisms) as a purely algebraic condition on 
their corresponding algebraic curvature tensors. Notice that the commuting of 
Killing tensors as endomorphisms on the tangent space is not equivalent to the 
commuting of their respective algebraic curvature tensors as endomorphisms on 
A^V^, as one may guess naively. 

Proposition 5.3 (Algebraic counterpart of the commutator). Let K and K he two 

Killing tensors on a non-fiat constant curvature manifold with algebraic curvature 
tensors R and R. Then the following statem,ents are equivalent 

[K, K]^0 (5.4a) 

I — \ — I — \ — I * 

bi 62 di d2 

- (5.4b) 



Q2 
C2 



9ijR^bia2b2^'' diC2C 







\bi\b2\di\d2\ giiR\^a^h2^ diC2d2 " ^ 



(5.4c) 



Proof. Using V^x'' 
V = da and w = dg: 



we write (1.6) in local coordinates for coordinate vectors 



The product of K and K, regarded as endomorphisms, is then given by 

Ra\h\a2h2^c\d\C2d2'^ XXX "V X V^X V'^X V^X 
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In [Schl2] we proved the following identity: 

As a consequence of the antisymmetry of algebraic curvature tensors in the last 
index pair, the term does not contribute when substituting this identity into 

the previous expression: 

K''^K\ =g''-'''Ra,b,a.MRc,d,c^d,x''^x''^x'^x'-V''x'''Vpx''^ . 
The commutator [K, K] is therefore given by 

[K,kU,^g''^''^Ra,b^a,b,Rc,d^c,d,x''^x''^x^^x''^W[^x''^Vp^x''^ 
and vanishes if and only if 

g''''''Ra,b,a,b,Rc,d,c,d,x'''x''^x''x'^v^'''w''^^ =Q (5.5) 
for all X E M and v,w £ T^M. That is, for all x,v,w E V with 

g{x,x) = l g{x,v) = g{x,w) ^0. (5.6) 

We can drop the restriction g{x,x) ~ 1 since M.M C V is open. We can also 
drop the restrictions 17(2;, v) = and g{x, w) = by decomposing arbitrary vectors 
v,w E V under the decomposition V = T^M (B^x. To see this, notice that (5.5) is 
trivially satisfied for w = a; or for w = x. Indeed, in this case the tensor 



Raibia2b2 RcidiC2d2 



(5.7) 



is implicitely symmetrised over five indices and Dirichlet's drawer principle tells 
us that this comprises a symmetrisation in one of the four antisymmetric index 
pairs. This means we can omit the restrictions (5.6) completely. In other words, 
[K,K] = is equivalent to (5.5) being satisfied for all x,v,w E V. 
Now notice that the tensor 

x''^x''^x'''x''^v^'''w'^^^ 

appearing in (5.5) is completely symmetric in the indices ai, 02, Ci and C2 and 
completely antisymmetric in the indices bi and d2. Applying (5.3) to this tensor 
therefore yields 



x''^x''^x''^x''^v^'''w'^^^ 



1 

2! 



d2 



1 



^ ' TT k^i P2 (^2 X X X V W 



a-i ^a2 ci C2 [bi d2] 



4! 



Ql 


02 \c1\c2\ 


a I 


02 Cl C2 


61 




61 




d2 




ch 





x'''x''^x''^x''^v^''^w'^^^ 



^ h ffli|g2|ci|c2 fai ai |a2 I ci I C2 | ^ai ^Ci ^C2 ^,[bi ^,^2] 

34560 y y 

(5.8) 



From the definition of the Young tableaux and their adjoints 

Ql 





02 Cl C2 1 






d2 





|ai |a2|ci |C2 



d2 



ai 


02 ci 1 <:2 


bi 









gl|"2 |Ci |C2 
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together with the properties 

|a2 I Ci I C2 I XXX — ^;x XXX 

and 





1 








ai 


bi 


= 3! 








da 



hi |ai |a2 I Ci |c2 — 5! hi |ai | 02 | Ci | C2 



we see that (5.8) simpHfies to 



ai 


02 ci 






d2_ 





+ 



2-5! 



\ 10368 m ' 34560 

When substituted into (5.5), this yields 



^SEEZH* ) x'''x''''x'''x''^v'''w'^^ 



ai 


02 Cl 1 C2 


h 




d2_ 





1 


h 


Ol a2|ci C2 


144 


d2_ 





As elements of the group algebra of the permutation group, the Young tableaux are 
linear combinations of permutations tt of the indices ai, a2, Ci, C2, 6i, (i2- For any 
such permutation tt we have 

R ,, , fi ., , 7.^(ii)T,^(a2)^'r(ci) 7r(c2) 7r(hi) 7r((i2) 
y 1-010102 02 ^ ''Ci aiC2 02 X X X X f 

= a''^'^^B ./ ^ *i \h R *i \,i w ^ *t,i ^x°-^x°"^x''^x''^v^^w'^'' 

y ^7r* (oi )7r* (01 Jtt* (02 J02 TT* (ci joiTT* (C2 )7r* (02 ) X X X X L/ uy , 

where tt* = tt^^ is the adjoint of tt. We can thus replace the Young tableaux in 
(5.9) acting on upper indices by its ajoint acting on lower indices: 



72 



Ol 


02 1 Ci C2 


61 




d2_ 





_J_ &i | oi | q2 | ci | c 2] 1 fc2<il D . R . j 

j y J''Oi6iO262-'"'Cl0lC2(l2 



= (5.10) 

Now notice that the second Young tableau involves a symmetrisation over the five 
indices 61, ai, 02, ci, C2 and that, as above, the symmetrisation of the tensor (5.7) 
in any five indices is zero. Hence the second term in (5.10) vanishes and we obtain 



Ol 


02 Cl C2 1 


61 




d2_ 





g''"''Ra,b,a^b,Rc,d,c,d, I X'^^X^^^X'^X'^V'^W''-' ^0. 



Ol^02^Ci C2,,^>l„,,d2 



Recall that \K, if ] = is equivalent to this condition being satisfied for all x, w G 
V . By polarising in x we get 



Ol 


02 Ci C2 


61 




d2 





g'"'''Raibia,b,Rc,dic,d,j x'^^y'^^Z^H'-v'^w"^ = 

for all X, y, z^t^v^w ^ V and hence 

rfi^di p p f\ 



ai^,a2 ^ci>C2„,bi d2 



01 


02 Ci C2 


61 




d2 
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This is the same as (5.4b) after appropriately renaming, lowering and rising indices. 
We have proven the equivalence (5.4a) <^ (5.4b). 

We now prove the equivalence (5.4b) -i^ (5.4c). Start from (5.4b) by expanding 
the Young tableau: 



a2_ \bi\b2\di\d2\ giiK 6^0262^ diC2d2 



0. 



In order to sum over all 4! permutations when carrying out the symmetrisation in 
the indices 61, &2, (^2, one can first take the sum over the 4 cyclic permutations 
of &i, 62, di, ^2, then fix the index bi and finally sum over all 3! permutations of the 
remaining 3 indices b2,di,d2'. 



\ b2 \ d1 \ d2 \gij (^R\^a,^b2^\iC2d2 + ^^b2a2diR\2C2h 
+ ^ dia2d2^'' b^C2b2 ^ £(202^1 ^"'b2C2di 







For a better readability we underlined each antisymmetrised index. Permuting 
the indices of the terms in the parenthesis under symmetrisation in &2 , c'l , '^2 and 
antisymmetrisation in &i , 02 , C2 we can gather the first and last as well as the second 
and third term: 



a2_ \ b2 \ di \ d2 \gtj[{R\^a,^b^ + ^^a^fci )-R''c;iC2d2 



" R\ia2d2 i^\^C2b2 + ^\2C2b^ ) ) — 



Using the symmetries of Rbia2b2i terms in the inner parentheses can be rewritten 
as 



pi E>1 

^ 0.26261 61620,2 



resulting in 



\b2\d1\d2\ gii 



0,26261 -^"'diC2(i2 ~^ ^ dia,2<i2 ^"'£36261) ~ ^ ■ 



As above, when carrying out the antisymmetrisation over &i,a2,C2, we can first 
sum over the three cyclic permutations of 61, a2, C2, then fix hi and finally sum over 
the two permutations of 02 , C2 . This results in 

I 62 I rfl I I gij (^'02 6261 ^''di C2 d2 + -^^rfl 02^2 -^^ 026261 



^ C262a2-^"'rfi6id2 ^ ^ diC2d2^\ib2a2 
^\ib2C2^^ dia2d2 + '^*<ii6i<i2 -^■'026202 



= 0. 



The symmetrisation of this in fei , 62 , di , 1^2 yields 



b2\di\d2\ gij[R\^j,^hiR 



1 f>i tdJ 

£2^2 ' dia2rf2 £26261 



+ ^'02620.2 ^''di6i(i2 + ^^diC2d2^\ib2a^ 



^ 616202-^ dia2(i2 



pi dJ 

(ii6i(i2^ O262C2 



= 0, 
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where we again underlined antisymmetrised indices. By the antisymmetry of alge- 
braic curvature tensors in the second index pair, all but the fourth and fifth term 
in the parenthesis vanish and we get 

I h I &2 | di I rf2 I ffij {R\ic^d2 ^'^bia^&a ^ ^^bic^b2^^ dia^d^ = ^ 

Due to the symmetrisation and antisymmetrisation we can permute the indices of 
the terms inside the parenthesis to find that 

fti I 62 | di | rf2 I 9r'jR\ia^b2 dtc^d2 ^ 0. 

Recall that this equation has been obtained from (5.4b) by a symmetrisation. This 
proves (5.4b) (5.4c). The converse follows easily by antisymmetrising (5.4c) in 
6i , 02 , C2 . This achieves the proof of the theorem. □ 

Remark 5.4. Obviously (5.4b) must be true for R = R. This fact is not evident, 
but has been proven in [Schl2j. 

For diagonal algebraic curvature tensors R and 7? on a four-dimensional vector 
space the condition (5.4c) can be written in the form 





/I 


Rijij 


Rijij \ 




det 


1 


Rjkjk 


Rjkjk = 


(5.11) 






Rkiki 


Rkiki / 





for all distinct i, j, k £ {0, 1, 2, 3}. 

Corollary 5.5. Stdckel lines'^ in the KS-variety correspond to Stdckel systems. 
More precisely, the preimage of a Stdckel line under the map from integrahle Kil- 
ling tensors with diagonal algebraic curvature tensor to the KS-variety is a Stdckel 
system. In particular, the Killing tensors in a Stdckel system have simultaneously 
diagonalisable algebraic curvature tensors. 

Proof. A Stackel line is the projective line through the (coUinear but not coinciding) 
points M, t(7r(M)) and ^(7r(M)) in the KS-variety for any non-singular KS-matrix 
M. We have to prove that corresponding Killing tensors Kq, Ki and K2 mutually 
commute. We first show that Ki and K2 commute. Let R and R be their diagonal 
algebraic curvature tensors. From their KS-matrices (4.8) and = — we 
can read off the values 

Wa =0 ta = Ua Wa = n\ t q = UpTL^ , 

where we have neglected the scalar curvature, which can be chosen appropriately. 
The diagonals of R and R are then given by the columns of the matrix 



/l 


Rqwi 


-Roioi^ 




(l 


Wi 


+ tl 


Wl 


+ ?1^ 
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+ni 


9 

nf 


+ n2n^ 
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-R0202 
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+ t2 


W2 


+ ?2 




1 


+n2 




+ n^ni 


1 


-R0303 


-^0303 




1 


W3 


+ t3 


W3 


+ h 




1 


+"3 


9 

ni 


+ nin2 


1 


^2323 


-R2323 




1 


Wi 


~ti 


Wl 
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-ni 


9 

nf 


- n2n3 


1 


^3131 


-R313I 




1 


W2 


-t2 


W2 






1 


-n2 


9 


- n-iUi 


V 


^1212 


-R1212 / 




V 


W3 


-h 


W3 


-k) 




V 


-713 


9 

ni 


- nin2j 



^ci. Definition 4.10 



a2 

C2 
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The four square matrices in (5.11) can be obtained from this matrix by discarding 
the rows (1,2,3), (1,5,6), (2,4,6) respectively (3,4,5) and it is not difficult to see 
that their determinants vanish. This shows that Ki and K2 commute. 

If t(7r(Af)) and u{t:{M)) are linearly independent, then AI is a linear combination 
of both and hence Kq is a linear combination of g, Ki and K2. In other words, Kq, 
Ki and K2 mutually commute. 

If l{'k{M)) and i/(7r(M)) are linearly dependent, then n = 7r(M) satisfies \ni\ — 
1*^21 — I "-si- For simplicity, suppose ?ii = 71,2 = 713 = 1. The other choices of the 
signs are analogous. By Lemma 4.7, M is of the form (4.6a) and thus has Wa — —ta- 
Consequently, the diagonal algebraic curvature tensors R and R of Kq respectively 
K2 are given by 
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^0101 


-Roioi^ 
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Wl 


+ tl 


Wl 
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2\ 
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^0202 


-R0202 
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+ t2 
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+ ?2 




1 
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^0303 


-Ro303 
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+ t3 


W3 


+ ^3 
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^2323 


-R2323 
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Wl 


-ii 


Wl 


-ii 
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2ti 
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^^3131 


-R313I 
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W2 


-t2 


W2 


-i2 
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2t2 





V 


^^1212 


-R1212 J 




V 


W3 


-h 


W3 


-k) 




V 


2<3 


V 



As above, they satisfy the commutation condition (5.11) and we conclude that Kq, 
Ki and if 2 commute. □ 

We can now prove theorem Theorem 1.10. Recall from Proposition 4.8 that 
each Stackel line contains a unique skew symmetric point i(n), where n e P^, and 
that each non-singular skew symmetric point t(n) determines a unique Stackel line, 
given by the two points b{n) ^ v[n). Moreover, the singular skew symmetric points 
are the four points n — (±1 : ±1 : ±1) € for which t(n) — v{n). Hence it suffices 
to show that the subspace spanned by t(n(i)) and i'{n{t)) has a well defined limit 
for t — > if n(0) is one of these and that this limit depends on ?i(0) G P^. We leave 
it to the reader to verify that the limit space is spanned by /,(n(0)) and /,(ri(0)). □ 

6. Geometric constructions of integrable Killing tensors 

In this section we present several geometric constructions of integrable Killing 
tensors and interpret each of them within our algebraic picture from the last section. 

6.1. Special Killing tensors. There are two ways in which integrable Killing ten- 
sors are related to geodesically equivalent metrics. The first is via special conformal 
Killing tensors. 

Definition 6.1. A special conformal Killing tensor on a Riemannian manifold is 
a symmetric tensor Lap satisfying the Sinjukov equation 

V^Lap — Xagpj + ^pgaj (6.1a) 

where 

A=iVtrL, (6.1b) 
as can be seen from contracting a and (3. 

Special conformal Killing tensors parametrise geodesically equivalent metrics in 
the following way [Sin66]. 
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Theorem 6.2. A metric g is geodesically equivalent to g if and only if the tensor 

( detg' 

is a special conformal Killing tensor. 

From the definition we see that special conformal Killing tensors on a manifold 
M form a vector space which is invariant under the isometry group of M . In other 
words, they define a representation of this group. The following lemma shows that 
this representation is isomorphic to a subrepresentation of the representation of the 
isomorphism group on Killing tensors. 

Lemma 6.3. If L is a special conformal Killing tensor, then K := L — (tv L)g is 
a Killing tensor. This defines an injective map from the space of special conformal 
Killing tensors to the space of Killing tensors, which is equivariant with respect to 
the action of the isometry group. 

Proof. It is straightforward to check the first statement. Taking the trace on both 
sides of K = L — {trL)g yields 

tri^ = (1 - n)trL 

and hence 

L = K-^g. (6.3) 

n — 1 

This shows injectivity. Equivariance is obvious. □ 

Definition 6.4. We define a special Killing tensor to be a Killing tensor of the 
form 

K = L-{tTL)g, 
where L is a special conformal Killing tensor. 

Proposition 6.5. A special conformal Killing tensor as well as the corresponding 
special Killing tensor have vanishing Nijenhuis torsion. In particular, special Killing 
tensors are integrable. 

Proof. Substituting the Sinjukov equation (6.1) into the expression (1.3) for K ^ L 
shows that the Nijenhuis torsion of L is zero. Substituting K = L — {tr L)g into 
(1.3) shows that this implies that the Nijenhuis torsion of K is also zero. □ 

For a non-flat constant curvature manifold M C 1^ we can identify the sub- 
representation of special Killing tensors inside the space of Killing tensors. Recall 
that the space of Killing tensors is isomorphic to the space of algebraic curva- 
ture tensors on V and an irreducible GL(y)-representation. Under the subgroup 
SO(y) C GL{V) this space decomposes into a Weyl component and a Ricci com- 
ponent. 

Definition 6.6. The Kulkarni-Nomizu product of two symmetric tensors h and k 
on a vector space V is the algebraic curvature tensor h & k, given by 

(/l ® k^ aibia2b2 • ^aia2^??ife2 ^0162^^)102 ^6102^0162 ^ ^6162^0102 ■ (^■^) 
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The Ricci component consists of all algebraic curvature tensors of the form 
where /i is a symmetric tensor. Therefore the dimension of the Ricci part is equal 
to N{N + l)/2 where N — dimT^. This is exactly the dimension of the space 
of special conformal Killing tensors on M [Sin66]. But there is no other SO(l^)- 
subrepresentation of the same dimension inside the space of algebraic curvature 
tensors. This shows the following. 

Lemma 6.7. A special Killing tensor on a non-flat constant curvature manifold is 
a Killing tensor whose algebraic curvature tensor R has a vanishing Weyl tensor. 
That is, R is of the form 

R^ h® g 

for some symmetric tensor h on V, where "®" denotes the Kulkarni-Nomizu pro- 
duct (6.4). 

Alternatively, this lemma can be checked directly by verifying that (6.3) satisfies 
(6.1) if K has an algebraic curvature tensor of the form R = h ® g. 

The following consequence is a special case of the so called "cone construction" 
[MMIO]: Any special conformal Killing tensor on a Riemannian manifold M can 
be extended to a covariantly constant symmetric tensor on the metric cone over M . 
In our case, where M C F is a constant curvature manifold, this cone is nothing 
but the embedding space V . That is why a special conformal Killing tensor is the 
restriction of a constant symmetric tensor on the ambient space. 

Proposition 6.8. Let R — h ® g be the algebraic curvature tensor of a special 
Killing tensor K = L — (tri)g on a non-flat constant curvature manifold M <ZV . 
Then the corresponding special conformal Killing tensor 

iiK , , 

L^K -g 6.5 

n ~ 1 

is the restriction of the (constant) symmetric tensor 

L = h-^g (6.6) 

from V to M , where "Tr" denotes the trace on V and "tr" the trace on the tangent 
space. In particular, if h is tracefree, then L is simply the restriction of h from V 
to M. 

Proof. The special Killing tensor K is given by substituting R = h ® g into (1.6) 
using (6.4). This yields 

Kx{v, w) = h{x, x)g{v, w) — h{x, v)g{x, w) — h{x, w)g{x, v) + h{v, w)g{x, x) 

for tangent vectors v,w Cz T^M at a point x G M, i.e. vectors x,v,w € V with 
g{x, x) = 1 and g{x, v) = g{x, w) — 0. This simplifies to 

Kx{v, w) — h{v, w) + h{x, x)g{v, w) . (6-7) 

Let ei, . . . , e„ be an orthonormal basis of T^M and complete it with cq := x to an 
orthonormal basis of V. Then 



tr/v = ^ Kx{ea,ea) = ^{h{ea,ea) + /i(eo, eo)g(e„, ea)) 

Q — 1 a— 1 

n 

— /i(ei, e^) + (n — l)/z(eo, e^) = Tv h -\- (n — l)h{x, x) 
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Recall that "tr" denotes the trace on T^M, whereas "Tr" denotes the trace on V . 
Substituting (6.7) and (6.8) into the right hand side of (6.5) now yields (6.6). □ 

For we can verify the integrability of a special Killing tensor on alge- 
braically, using the Theorems 1.4 and 1.5: 

Lemma 6.9. (i) R — h <& g is diagonalisahle. More "precisely, R is diagonal 
in a basis where h is diagonal. In this basis the diagonal elements of R are 
given in terms of the diagonal elements hi := ha by 

Rijtj = /^^ + h-j j ■ 

(ii) If divaV — 4, the KS-matrix (1.9) of R — h® g is antisymmetric with 

Aa =0 ta = ho + ha ■ 

In particular it has determinant zero. 

Corollary 6.10. The space of special Killing tensors on corresponds to the 
projective space of antisymmetric matrices inside the KS-variety. 

Proposition 4.8 and Corollary 5.5 now give a unique representative in each Stackel 
system. 

Corollary 6.11. Every Stackel system on contains a special Killing tensor 
which is unique up to multiplication with constants and addition of multiples of the 
metric. 

6.2. Killing tensors of Benenti type. The second way in which integrable Kil- 
ling tensors arise from geodesically equivalent metrics is described by the following 
theorem [LC04, Pai97, Bcn92, MT98]. 

Theorem 6.12. If a metric g on a (pseudo-)Riemannian manifold is geodesically 
equivalent to g, then 



\detg 



is an integrable Killing tensor for g. 



Definition 6.13. We will call a Killing tensor of the form (6.9) a Benenti- Killing 
tensor. 

Recall that the standard metric on the unit sphere S*" C in a Euclidean vector 
space (y, g) is the restriction of the scalar product g from F to 5" , which we denoted 
by g as well. Its geodesies are the great circles on S*". Consider the map 

/: S'' 

for A € G1j{V). Since A takes hyperplanes to hyperplanes, / takes great circles to 
great circles. Hence the puUback g of g under / is geodesically equivalent to g by 
definition. The Killing tensor obtained from applying (6.9) to g is given by [Schl2] 

g{Ax,Ax)g{Av,Aw) - g{Ax,Av)g{Ax,Aw) ^ rr. 

(detA)^ 
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Comparing with (1.6), we see that the algebraic curvature tensor of this KiUing 
tensor is proportional to h ® h, where "®" denotes the Kulkarni-Noniizu product 
(6.4) and 

h{v, w) := g{Av, Aw) = g{v, A^-Aw) . 

Moreover, we can replace the scalar product h by any - not necessarily positive 
definite - symmetric tensor. This motivates the following definition. 

Definition 6.14. We will say a Killing tensor on a constant curvature manifold of 
Benenti type if its algebraic curvature tensor is proportional to 

R = h&h 

for some symmetric tensor h on V . 

A Killing tensor of Benenti type with positive definite /i is a Benenti-Killing 
tensor and thus integrable. With the aid of the algebraic integrability conditions 
one can verify that any Killing tensor of Benenti type is integrable [Schl2]. On 
we can check this directly, using the Theorems 1.4 and 1.5: 

Lemma 6.15. (i) For a symmetric tensor h on a vector space V , the algebraic 
curvature tensor R — h(S)his diagonalisahle. More precisely, R is diagonal 
in a basis where h is diagonal. In this basis the diagonal elements of R are 
given in terms of the diagonal elements hi := hu by 

R,j,j^2h,hj i^j. (6.10) 

(ii) If dim V = 4, the KS-matrix (1.9) of R~ hQSih is given by 

Aq = (/lo - ha){hp - h^) ta = hoha - h/sh^y 

and has zero determinant and trace. 

6.3. Benenti systems. Special Killing tensors and Benenti-Killing tensors only 
constitute particular cases of a more general construction of integrable Killing ten- 
sors out of special conformal Killing tensors. This construction appears in the 
literature in different guises and under a number of different names: "Newtonian 
systems of quasi-Lagrangian type" [RWML99] , "systems admitting special confor- 
mal Killing tensors" [CSTOO, CraOSa], "cofactor systems" [LunOl, RW09], "bi-quasi- 
Hamiltonian systems" [CS02, CraOSb] or "L-systems" [Ben05]. Here we follow Bolsi- 
nov & Matveev, who call them "Benenti systems" [BM03], because their formulation 
is directly reflected in our algebraic description for S^. 

Recall that the adjugate matrix Adj M of a matrix AI is the transpose of its 
cofactor matrix, which is equal to (det M)M~^ if M is invertible. Given a metric, 
we can identify endomorphisms and bilinear forms. This allows us to extend this 
definition to bilinear forms. As for endomorphisms, we denote the adjugate of a 
bilinear form L by Adj L, which is again a biliner form. 

Theorem 6.16 ([BM03]). For a special conformal Killing tensor L the family 

K{\) ^ Adi{L- Xg) AeM 

is a family of mutually commuting integrable Killing tensors. In particular, for 
X = we recover the Benenti-Killing tensor (6.9) from (6.2). 

Definition 6.17. The family K{X) = Adj(L — Xg) of mutually commuting Killing 
tensors is called the Benenti system of the special conformal Killing tensor L. 
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In order to state the algebraic counterpart of the above theorem for constant 
curvature manifolds, we first need the following lemma. 

Lemma 6.18. Let h be a symmetric tensor on a vector space V of dimension n + 1. 
Then the algebraic curvature tensor 

(Adjfe)®(Adj/.) 

deth ^ ' ' 

is well defined and homogeneous in h of degree n — 1. 

Proof. Chose a basis in which h is diagonal with diagonal elements hi := ha. Then 
Adj h is also diagonal and has diagonal elements ho ■ ■ ■ hi-ihi^i ■ ■ ■ hn- By (6.10) 
the algebraic curvature tensor (6.11) is diagonal and given by 

ho ■ ■ ■ hi^ihij^i ■ ■ ■ hn ■ ho ■ ■ ■ ■ ■ • 



— 2/io • • • ft-i-i/ii+i • • • hj^ihj^i ■ ■ ■ hn 
for i ^ j. This is well defined and homogeneous in h of degree n — 1. □ 

Theorem 6.19 (Algebraic representation of Benenti systems). Let L be a special 
conformal Killing tensor on a non-fiat constant curvature manifold M C V , i. e. 
the restriction of a (constant) symmetric tensor L from V to M. Then the algebraic 
curvature tensors of the Benenti system 

are given by 

^ Adj(£- Aff)® Adj(£- A.g) 
det(L-Ag) 

In particular, a member of this Benenti system is of Benenti type unless A is an 
eigenvalue of L and it combes from a geodesically equivalent metric if —A is suffi- 
ciently large. 

Proof. For simplicity of notation, set h := L — Xg and A := Adj h and let K be the 
Killing tensor corresponding to the algebraic curvature tensor R{X). Substituting 
R{X) = {A® A) /deth into (1.6) using (6.4) yields 

, , A(x,x)A(v,w) — A(x,v)A(x,w) 

K(V,W) = ; ; 

^ deth 

for v,w G T^M. Our aim is to show that K — Adj(L — Xg). We can extend the 
bilinear form K from T^M = C 1^ to a bilinear form on 1/ by allowing any 
v,w €V in this expression. In coordinates, this reads 

j> Akix^x^Aij — AikX^ Ajix^ 



det h 



Notice that K{v, w) ~ ior v — x or w ^ x. This means that K is the extension 
of if by under the orthogonal decomposition V = x^ (B Mx. Therefore we have 

K • (L - Ag) = KU. • (L - Ag)U. = KU. • hj,. = (K • h)U. , (6.12) 
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where the boldface letters indicate that we now consider endomorphisms. By the 
definition of A we have A^jh\^ = (det h)g'-„^, so that 



{Akix''x'A'j - A\x^Ajix^)h\ 



= ^kix^x^ g^jyi — A^ j^x^ gmix^ ■ 

Recall that v S T^M means g{x,v) = 0. Consequently the second term vanishes 
on T^M and we get 

(K-h)U^ =^(x,x)g-(dethU^)g = det(L-Ag)g. (6.13) 

Here the second equality follows from the definition of the adjugate as the transpose 
of the cofactor matrix. Combining (6.12) with (6.13) we obtain 

K.(L-Ag)=dct(L-Ag)g 

for all A e M and hence K — Adj(L — Xg). This accomplishes the proof. □ 

The adjugate of an (ri x rt)-matrix M is a polynomial in M of degree n—l, given 

by 

Xm(A) - Xa/(0) 
A 

where xm is the characteristic polynomial of M. Therefore a Benenti system is a 
polynomial in A of degree n — 1: 

n-l 

if(A)=Adj(L-Ag) = ^if,A^ 

Up to signs, the leading coefficient of this polynomial is the metric, the next coeffi- 
cient is the special Killing tensor (6.2) and the constant term is the Benenti-Killing 
tensor (6.9): 

Kn-i = lim -T^K{X) = lim Adj(iL - g) = (-l)"-ig 

A— !-oo A A— i-oo 
d 



^"-^ dX 



A=0 "'^ 



Adj(Ai-.g)-(-l)"-i(L-(trL)5) 

A=0 



Ko = KiO) = Adj i . 



Note that since all the K{X) commute, the n coefficients Kj span a vector space 
of mutually commuting integrable Killing tensors. If the n coefficients are linearly 
independent, then this space is a Stackel system and the Kj form a basis therein. 
However, it is in general an open question under which conditions this will be the 
case. For S'^ we can derive the answer directly from our algebraic description. 

Lemma 6.20. For a symmetric tensor h on a four dimensional vector space we 
have 

(Adj /i) © (Adj /i) ,s ,^ , 

— ^—L — ^ — ^=*(/i©/i*, (6.14) 
det n 

where * is the Hodge star operator. 

Remark 6.21. Regarding (2.2), the conjugation of an algebraic curvature tensor 
with the Hodge star is simply a sign change in the trace free Ricci part: 

R = W + T + S *R* = W T + S 
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Proof. Choose a basis in which h is diagonal with diagonal elements hi := ha. For 
simplicity of notation we assume that det h ^ 0. Then hi ^ and Adj h is diagonal 
with diagonal elements Aeih/hi. Substituting Adj h for h in (6.10) and comparing 
to (6.10) shows that the diagonal algebraic curvature tensors 

R = (Adj h) ® (Adj h) R = h®h 

are related by 

~ ^dethdeth ,x ^0/11^2^3 , , /, , \ 

Rijij = ^—^ ^ = 2(det/i) — j-j- — = [deth)2hkhi = [dethjRkm , 

for {i,j, k, 1} = {0, 1, 2, 3}. Now the identity (6.14) follows from (2.13) and the fact 
that R and *R* differ in the sign of the trace free Ricci part. □ 

Proposition 6.22. Let L be a special conformal Killing tensor on C V , i. e. 
the restriction of a (constant) symmetric tensor L from V to , and consider the 
corresponding Benenti system 

K{X) = Adj(L - \g) = K2\^ + KiX + Kq. 

Then the coefficients Kq, Ki and K2 span a Stackel system if and only if L has no 
triple eigenvalue. 

Proof. By Theorem 6.19 and (6.14) the algebraic curvature tensor of K{\) is the 
polynomial 

R{\) = *((L - \g) @{L- A5))* 

= ^^9 ® .9 + 2AL ® g + *{L ® L)* 

=: A^i?2 + Ai?i + i?o ■ 

The coefficients Rq, Ri and i?2 are the algebraic curvature tensors of Kq, Ki and 
K2 respectively. We can assume them to be diagonal. Since K2 is minus the metric, 
Kq, Ki and K2 are linearly independent if and only if Rq and Ri define two distinct 
points in the KS-variety. 

Choose a basis in which L is diagonal with the eigenvalues A^ on the diagonal. 
Without loss of generality we can assume that L is trace free. Motivated by (2.12), 
we can thus parametrise L by a vector n S as 

^ + n, + ^ ri^-n,-n, ^^^^^^^ 

where 

n„ = Ao + A„ = -(A^ + A.^). (6.15b) 

Then L has a triple eigenvalue if and only if \ni\ = \n2\ = \n:i\. This condition is 
equivalent to the condition that the KS-matrices i'{n) and t(n) in (4.8) are pro- 
portional, i.e. define the same point in the KS-variety. Our proof will therefore be 
finished if we show that Rq and Ri define the points v{n) respectively t(n) in the 
KS-variety. 

That the algebraic curvature tensor Ri ~ 2L & g defines the point t(n) is the 
statement of Lemma 6.9. And since Rq — *{L ® L)* differs from L ® L in the sign 
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of the trace free Ricci part, Lemma 6.15 implies that for Rq the KS-matrix (1.9) is 
given by 

Aa = (Ao - A„)(A0 - A^,) = - 
= -(AoAa - A/jA-y) ^npn^. 
Now compare this to v{n) in (4.8) to finish the proof. □ 

The Benenti system of a special conformal Killing tensors L defines a map M — > 
K.{S^) given by A i— >■ K{X) — Adj(i — \g). In the proof above we have actually 
shown that this defines a map from = M U {oo} to the KS-variety which maps 
I— 7- vin) and cxi i— >■ L{n), where n £ is determined by L. If the points i(n) and 
i'{n) coincide, the image of this map is a single point. If they are distinct, the image 
is - by definition - the Stackel line through both points. In this case the image of 
the Benenti system is this Stackel line without the point i{n). 

We can now prove Theorem 1.8. It suffices to show that the image of diagonal 
algebraic curvature tensors of the form h® h is dense in the KS-variety. By Theo- 
rem 6.19 the algebraic curvature tensors of this form are dense in a Benenti system, 
so it suffices to show that the image of Benenti systems in the KS-variety is dense. 
But this follows from the preceding paragraph and part (ii) of Proposition 4.8. □ 

6.4. Extension of integrable Killing tensors. If we fix an orthogonal decom- 
position IR"i+"2+2 _ K"2+i^ then two algebraic curvature tensors i?i on 
and i?2 on R"2+i define an algebraic curvature tensor i?i ® i?2 on M"i+"2+2^ 
Correspondingly, two Killing tensors Ki on S^'^ and K2 on 5*"^ define a Killing 
tensor Ki ® K2 on 5'"i+"2+i^ jj^ particular, we can extend any Killing tensor K on 
5'""^ by zero to a Killing tensor K ®Q on S*". 

From the algebraic integrability conditions we see that Ki K2 is integrable if 
Ki and K2 are both integrable. This direct sum operation defines an embedding 

JC{S"^) X ^(5"^) ^(5"i+"^+i) 

where IC{S") denotes the space of integrable Killing tensors on S"". In particular 
we have an embedding 

/C(S'"^i) -^/C(S'"). 

Definition 6.23. We will call (integrable) Killing tensors arising in the way de- 
scribed above extensions of Killing tensors from (S'"i,S'"^) respectively S"~^. 

As an example, let us determine which integrable Killing tensors on are 
extensions from respectively {S^,S^). Note that in dimension one and two all 
Killing tensors are integrable. This follows directly from the integrability conditions. 

An algebraic curvature tensor R on M"' is an extension by zero of an algebraic 
curvature tensor on M.^ under the inclusion {xi, X2, X3) 1— > (0, xi, X2, X3) if and only 
if Raibia2b2 zero whenever one of the indices is zero. In particular, the diagonal 
algebraic curvature tensor (1.8) on is an extension from 5^ if and only if Wa + 

— —ta- In this case the KS-matrix (1.9) is equal to (4.6a) with ni = n2 — n^. 
Therefore the extensions of Killing tensors from S'^ correspond to the isokernel 
planes in the KS-variety. In particular, the face centers of the octahedron in the KS- 
variety correspond to the extension of the metric on S"^ under the four embeddings 
5^ C given by intersecting 5''^ C M"* with each of the four coordinate hyperplanes. 

An algebraic curvature tensor R on M.^ is a sum under the decomposition of 
{xo,Xi,X2,X3) into (xqjXi) and {x2,X3) if and only if the components i?aibia2&2 
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are zero unless all four indices are either in {0, 1} or in {2, 3}. In particular, a 
diagonal algebraic curvature tensor (1.8) is an extension from (S^, S^) if and only 
ii W2 + j2 ~ t2 — ws + j2 ~ t-s ~ Q. In this case the KS-matrix (1.9) is equal to 



where V+i and V-i are rank one singular points in the KS-variety, c. f. (4.4). There- 
fore the extensions of Killing tensors from (5^, S^) correspond to the projective lines 
joining opposite vertices of the octahedron in the KS-variety. In particular, each 
vertex of the octahedron corresponds to an extension of the metrics on two ortho- 
gonal copies of in S^, given by intersecting C with a pair of orthogonal 
coordinate planes. 



We eventually demonstrate how separation coordinates on and their classifica- 
tion arise naturally and in a purely algebraic way from our description of integrable 
Killing tensors. 

Kalnins and Miller proved that separation coordinates on S" are always ortho- 
gonal [KM86]. From the work of Eisenhart we know that orthogonal separation 
coordinates are in bijective correspondence with Stackel systems [Eis34]. And in 
the preceding sections we have proven that on every Stackel system contains 
a special Killing tensor which is essentially unique. By definition, special Killing 
tensors are in bijective correspondence with special conformal Killing tensors. But 
Crampin has shown that the eigenvalues of a special conformal Killing tensor are 
constant on the coordinate hypersurfaces of the corresponding separation coordi- 
nates [Cra03a]. Hence, if these eigenvalues are simple and non-constant, they can 
be used as separation coordinates. This reduces the classification of separation co- 
ordinates on to a computation of the eigenvalues of special conformal Killing 
tensors. The computation is considerably simplified by the fact that a special con- 
formal Killing tensor on 5" C V is the restriction of a (constant) symmetric tensor 



7.1. Eigenvalues of special conformal Killing tensors. Let L be a special 
conformal Killing tensor on a non-flat constant curvature manifold M C V, i. e. 
the restriction from V to M of a (constant) symmetric tensor L. From now on 
we consider L and L as endomorphisms. We want to compute the eigenvalues 
Ai(a;) ^ . . . ^ A„(a;) of L: T^M — > T^M for a fixed point x G M and relate them 
to the (constant) eigenvalues Aq ^ . . . ^ An of L: V ^ V. 

Let P : V V he the orthogonal projection from V = T^M © to T^M. 
Then the restriction of PLP to T^M is L and the restriction of PLP to x is 
zero. Therefore x is an eigenvector of PLP with eigenvalue and the remaining 
eigenvectors and eigenvalues are those of L. This means the eigenvalue equation 
Lv — \v is equivalent to PLPv = Xv for u _L x. By the definition of P we have 
Pv — V and PLv = Lv — g{Lv,x)x. Hence we seek common solutions to the two 




7. Separation coordinates 



on V. 



equations 



(L — X)v = g{Lv, x)x 



g{x,v) = 0. 
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If A is not an eigenvalue of L, then {L — A) is invertible and v is proportional to 
(L — \)^^x. The condition v Jl x then yields the equation 

g{CL-\)-^x,x)^Q 

for A. In an eigenbasis of L this equation reads 

k=0 

That is, the zeroes of the function q{X) are eigenvalues of L. As depicted in Figure 2, 
q{X) goes to for A -> ±oo, has poles at the eigenvalues Ak of L unless a;^ = and 
is nionotonely increasing in between. 




Figure 2. Eigenvalues of a special confornial Killing tensor on 5". 



Let us first suppose that L has only simple eigenvalues. If Xk ^ for all k = 

0. . . . , n, then q has exactly n + I poles at Aq, . . . , A„ and n zeroes in between. If 
Xk = 0, then q has no pole at Aj.. But Xk = means that the k-th eigenvector of 
L is orthogonal to x and hence also an eigenvector of L with eigenvalue A^. We 
conclude that the eigenvectors Vk of L which are orthogonal to x are eigenvectors 
of L with the same eigenvalue Afe, and that the remaining eigenvalues of L are the 
zeroes A^ of the function q{X) with corresponding eigenvector = {L — \k)~^x. 
In particular, the eigenvalues of L and L intertwine: 

Ao < Ai(x) < Ai s$ X2{x) sc • • • sc A„(x) s$ A„ . (7.2) 

Note that the common denominator of the fractions in the sum (7.1) is the charac- 
teristic polynomial Xli^) of L. Therefore the characteristic polynomial of L is 

Xl(A)=-z(A)x£(A), (7.3) 

1. e. the function q is the quotient of the characteristic polynomials of L and L. If 
we write this quotient in the form 

n 

n 2 n (Afe — A) 

fc=o Afe A _ 

fe=0 

we see that it is the defining equation for elliptic coordinates Ai(a;), . . . , A„(a::) on 
5" with parameters Aq < Ai < • • • < A„, c. f. [WRW03]. 
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If L has a simple eigenvalue, then the corresponding eigenspace is one dimen- 
sional and has only a trivial intersection with a generic tangent space to Af , because 
Af C 1^ is of codimension one. If L has a double eigenvalue, then the corresponding 
eigenspace is a plane and has a one-dimensional intersection with a generic tangent 
space to A/. This defines a field of common eigenvectors for L and L on M . The 
corresponding eigenvalue is the same for L and L and hence constant on M. This 
also follows from (7.2). In the same way an eigenvalue of L with multiplicity m-\-l 
defines an m-dimensional eigenspace distribution of L with constant eigenvalue. 

Proposition 7.1. Let L he a special conformal Killing tensor on a constant curva- 
ture manifold M C V , i. e. the restriction of a constant symmetric trace free tensor 
L on V . Then: 

(i) If L has only simple eigenvalues, then the eigenvalues of L define elliptic 
coordinates on M whose parameters are the eigenvalues of L. 

(ii) Any multiple eigenvalue of L is a constant eigenvalue for L. 

(Hi) L has simple eigenvalues if and only if L has at most double eigenvalues. 

7.2. Killing tensors and separation coordinates on S*^. Before we turn to 
, let us briefly discuss integrable Killing tensors and separation coordinates on 
5^. The reason is that their extensions will appear on . First note that in 
dimension two every Killing tensor is integrable. So we do not have to deal with 
the integrability conditions. Note also that algebraic curvature tensors on have 
a zero Weyl part and are thus determined by their Ricci tensor T alone. This means 
that Killing tensors on 5*^ C M'^ correspond to symmetric tensors T on M"^, where 
the metric corresponds to the identity. 

Stackel systems on S"^ are of dimension two and contain the metric. That implies 
that each Stackel system on corresponds to a plane of Ricci tensors spanned by 
the identity and some trace free symmetric tensor, the latter being unique up to 
multiples. Moreover, the isometry group acts by conjugation on the Ricci tensor. 
Accordingly, the classification of separation coordinates on S"^ modulo isometrics is 
equivalent to the classification of trace free symmetric tensors T on under the 
orthogonal group 0(3). 

In particular, every Killing tensor on S"^ is special and we can make use of the 
results in the previous section. Depending on the multiplicities of the eigenvalues 
^1 ^ ^2 ^ ^3 of the trace free Ricci tensor T of a Killing tensor on S'^ we consider 
the following two cases. 

Elliptic coordinates. If ti < t2 < t^, then the restriction L of T from to C K'^ 
has two different eigenvalues, given by the zeros of its characteristic polynomial 
(7.3): 

Xl(A) = {t2 - X){t3 - X)xl + (t3 - A)(<i - X)xl + {ti - X){t2 - X)xl = . 

Using that xf + X2 + x'^ — I and ti + <2 + ^3 = and writing the endomorphisms 
T and Adj T as forms, this quadratic equation reads 

+ T{x, x)X + (Adj T){x, x)^0. 

The two solutions Ai(a;) and X2{x) are real and satify ti < Xi{x) < t2 < X2{x) < t^. 
As mentioned above, they define elliptic coordinates on with parameters ti < 
t2 < h. 
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type 


coordinates 


induced from 


[Eis34] 


[KM86] 


(0123) 


elliptic 




V 


(1) 


(01(23)) 


oblate Lame rotational 




I 


(2a) 


(0(12)3) 


prolate Lame rotational 




I 


(2b) 


((01)(23)) 


cylindrical 




III 


(5) 


(0(123)) 


Lame subgroup reduction 


S'2 elliptic 


IV 


(3) 


(0(1(23))) 


spherical 


spherical 


II 


(4) 



Table 1 . Classification of separation coordinates on 



Spherical coordinates. If ti < t2 — t^, the characteristic polynomial of L reduces to 

Xl{X) = (t2 - A) [{t2 - X)xl + ih - X){xl + xD] = . 
One of the eigenvalues only depends on xi and the other is constant: 

Xlix) ^ t2xl + ti{l ~ xl) X2ix)=t2 

The curves of constant Xi{x) are circles of constant latitude and the meridians 
complete them to an orthogonal coordinate system - the spherical coordinates. 

The case ti — t2 < is analogue and yields spherical coordinates as well. Note 
that even so, both coordinate systems are not isometric. This is because they differ 
not only by an isometry but also by an exchange of the two coordinates. 

We can use the action of the isometry group 0(3) to diagonalise the Ricci tensor 
and thereby the algebraic curvature tensor. Hence the projective variety of inte- 
grable Killing tensors on with diagonal algebraic curvature tensor is isomorphic 
to the projective plane P^, where the point (1 : 1 : 1) S corresponds to the metric 
on S^. Moreover, Stackel systems on correspond to projective lines through this 
point, since they always contain the metric. The residual isometry group action on 
Killing tensors with diagonal algebraic curvature tensor is the natural action of S3 
on P^ by permutations. 

7.3. The classification. We now recover the well known classification of separa- 
tion coordinates on S^ [Eis34, KM86]. In particular, this will explain the graphical 
procedure given in [KM86] in terms of the classification of trace free symmetric 
tensors under the orthogonal group. The results are summarised in Table 1. 

As before, we will denote by Aq ^ Ai ^ A2 ^ A3 the ordered eigenvalues of the 
symmetric tensor L associated to a Stackel system and parametrised by a vector 
(ni,n2,n3) S via (6.15a). Recall that we have the freedom to add a multiple of 
the metric to a Killing tensor or to multiply it with a constant. Both operations 
only affect the eigenvalues, not the eigenspaces, and therefore result in different 
parametrisations of the same coordinate system. Whereas Kalnins & Miller use 
this freedom to set Aq = and Ai = 1, here we use it to assume that L is trace free 
and then to consider (Ag : Ai : A2 : A3) respectively {ni : 77.2 : n^) as a point in P^. 

In particular, changing the sign of L reverses the order of the eigenvalues. This 
leaves us with the following five different alternatives for the multiplicities of the 
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ordered set of eigenvalues of L. For simplicity we write for "A^" and denote 
multiple eigenvalues by parenthesising them: 

(0123) (01(23)) (0(12)3) ((01)(23)) (0(123)). 

(0123) Elliptic coordinates. We have already seen that if the eigenvalues of L are 
simple, then the eigenvalues of L are also simple and define elliptic coordinates on 
5'" given implicitely by (7.4). 

(01(23)) Oblate Lame rotational coordinates. If L has a double eigenvalue Afc_i = 
Afc , then the corresponding eigenvalue of L is constant and can not be used as a 
coordinate. Nevertheless, we can derive the missing coordinate from the fact that 
the coordinate hypersurfaces are orthogonal to the eigenspaces of L. Let E be an 
eigenplane of L and denote by E-^ its orthogonal complement in V. Recall that 
the corresponding eigenspace of L at a point a; G 5"" is given by the intersection of 
E with the tangent space T^S"^ . Its orthogonal complement in T^S"^ is the inter- 
section of T^S"^ with the hyperplane that contains both x and E^ . The coordinate 
hypersurfaces are therefore the intersections of hyperplanes containing E^ with 5" 
and can be parametrised by the polar angle in E. Hence this angle can be taken as 
a coordinate. 

By (7.2), L has single eigenvalues given by (7.4) if and only if L has at most dou- 
ble eigenvalues and every double eigenvalue of L determines a constant eigenvalue 
for L. If we replace every constant eigenvalue by the polar angle in the correspond- 
ing eigenplane of L, we obtain a complete coordinate system. Concretely, in the 
case Aq < Ai < A2 = A3 the separation coordinates are defined by the solutions 
Ai(a;) and \2{x) of (7.4) together with the polar angle in the (x2, a;3)-plane. 

For Aq < Ai < A2 = A3 we have ni < 712 = 713. The corresponding Stackel line 
thus goes through the two points in the KS- variety given by the KS-matrices (4.8): 

771712 








-772 


n2 \ 


/ 


— nirt2 




"-2 





-Ul 




77l772 772 ~ ""1 


V 


-772 


771 


/ 


V 


-77,1772 772 



/ 





-772 


n2 




"-2 





-771 




-772 


771 







Subtracting 77i times the first matrix from the second, we see that the same Stackel 
line is also determined by the two points 



(0(12)3) Prolate Lame rotational coordinates. The case Aq < Ai = A2 < A3 is 
analogue and yields Lame rotational coordinates which differ from the above by an 
isometry and an exchange of two of the coordinates. 

((01)(23)) Cylindrical coordinates. In the case Aq = Ai < A2 = A3 Equation (7.4) 
has only one non-constant solution, namely A2(a::): 

Ai=Ai < \2{x)^ Ki{xl + xl) + K3{xl + x\) ^ A3 = A3 

The separation coordinates are obtained by replacing the constant eigenvalues Ai 
and A3 by the polar angle in the (xq^Xi) respectively (a:;2, a;3)-plane. 
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For Aq = Ai < A2 = A3 we have ni < n-z = = and the corresponding Stackel 
Hne goes through the two points in the KS- variety given by the KS-matrices (4.8): 



i.e. opposed vertices of the octahedron in the KS- variety. We have seen that this 
descibes extensions from {S^, S^) to S^. 

(0(123)) Coordinates extended from S"^ C . In the remaining case Aq < Ai = 
A2 = A3. Equation (7.4) has the solutions 



We can take the non-constant eigenvalue Ai(a;) as one separation coordinate on 
. However, since the other two eigenvalues of L are equal, we can not recover 

all separation coordinates from the special conformal Killing tensor L in this case. 

We will have to consider another Killing tensor in the Stackel system instead - one 

with simple eigenvalues. 

Before we do so, let us collect some facts from previous sections for the case 

Ao < Ai = A2 = A3. 

Proposition 7.2. Consider the following objects associated to a system of separa- 
tion coordinates on : 

• the Stackel system 

• a special conformal Killing tensor L whose eigenvalues are constant on the 
coordinate hypersurfaces 

• the Benenti system K{X) = Adj(L — \g) of L 

• the special Killing tensor K = L — {tv L)g in the Stackel system 

• the Stackel line in the KS-variety 

• the trace free symmetric tensor L on V which restricts to L 

• the vector n = (ni, 712,^3) € parametrising L via (6.15a) 
Then the following statements are equivalent: 

(i) The Stackel system consists of extensions of Killing tensors from S'^ C 
(up to the addition of a multiple of the metric). 

(ii) L has a multiple eigenvalue (which is constant). 
(Hi) if (A) does not span the Stackel system. 

(iv) K corresponds to a singular skew symmetric point in the KS-variety. 
(v) The Stackel line lies on an isokernel plane in the KS-variety. 
(vi) L has a triple eigenvalue, 
(vii) \ni \ = |n2| = |n3| 

Proof. The equivalence of (vi) and (vii) is a direct consequence of (6.15a). The 
equivalence to the other statements follows from (i) Section 6.4 (ii) Proposition 7.1 
(iii) Proposition 6.22 (iv) Proposition 4.4 and (v) Lemma 4.7. □ 




Equivalently, this Stackel line goes through the two points 




Ai(a;) = AiXg + Aoixj + xj + xj) 



A2 — A3 — Ai . 
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To make the extension of Killing tensors from S"^ C more explicit, consider 
a Killing tensor on S*^ with algebraic curvature tensor R and Ricci tensor T and 
choose a basis in which T is diagonal with diagonal elements ti , t2 and . Then R 
is also diagonal and, if T is trace free, given by 

-R1212 =^1+^2 = ^2323 = ^2+^3 = ^^1 -^3131 = ^3 + ^1 = ^^2 • 

In Section 6.4 we saw that this Killing tensor extends to a Killing tensor on 
corrsponding to the point in the KS-variety that is defined by the KS-matrix 




Vt-.^ h-ts -h \ = hV+i + t2V+2 + t3V+ 



If T is not trace free, we obtain the same result. This shows that the eigenvalues 
of the Ricci tensor of a Killing tensor on S'^ induce barycentric coordinates on the 
isokernel planes in the KS-variety. 

In Section 7.2 we have already computed the eigenvalues A2(a::) and Xsix) of a 
Killing tensor on S*^, parametrised by the eigenvalues of its Ricci tensor T. We see 
that if one perturbs the special Killing tensor on S'^ within a Stackel system, the 
constant double eigenvalue A2 = A3 splits up into two simple eigenvalues A2 (x) and 
A3(x), which come from separation coordinates on under an embedding S'^ C S'^. 
On the other hand, the eigenvalue Ai(a::) remains unchanged. Together, they define 
separation coordinates on S^. We denote this by writing the multiplicities of L 
as (0(123)) respectively (0(1(23))), depending on whether the multiplicities of the 
Ricci tensor T are (123) or (1(23)). This yields the remaining two systems of 
separation coordinates: 

(0(1(23))) Spherical coordinates. These are extensions of spherical coordinates on 
5^ and correspond to Stackel lines through a face center and a vertex of the octa- 
hedron in the KS-variety. 

(0(123)) Lame subgroup reduction. These are extensions of elliptic coordinates on 
and correspond to Stackel lines through a face center, but not a vertex of the 
octahedron in the KS-variety. 

7.4. The classification space. The non-generic separation coordinates are char- 
acterised by rii — ±n2, n2 = ±^13 and — ±ni. This defines six projective lines 
that divide the projective plane of skew symmetric matrices in the KS-variety into 
twelve congruent triangles. The four singular skew symmetric points are given by 
I rill — 1 1 = Irisl and each of the triangles has two of them as vertices. On the 
blow-up in these points, the triangles become pentagons. Any of these pentagons 
constitutes a fundamental domain for the S'4-action. Regarding our notation for 
the multiplicities of L, we realise that the vertices are labeled by the different 
possibilities to parenthesise the ordered set of eigenvalues of L with pairs of cor- 
rectly matching parentheses and that the edges correspond to a single application 
of the associativity rule. This naturally identifies the classification space for separa- 
tion coordinates on with the associahedron K4, also known as Stasheff polytope 
[Sta63, Lee89, Sta98]. It also shows that the space of Stackel lines in the KS-variety 
is naturally isomorphic to the permutoassociahedron KP4 [Kap93]. 
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(((01)2)3) 
spherical 



Lame 
subgroup 



(0((12)3)) ((0(12))3) 
spherical prolate spherical 

Lame 

Figure 3. Identification of the moduli space of separation coor- 
dinates on with the associahedron K4. 
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